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We survey 129 lepton number violating effective operators, consistent with the minimal Standard 
Model gauge group and particle content, of mass dimension up to and including eleven. Upon 
requiring that each one radiatively generates the observed neutrino masses, we extract an associated 
characteristic cutoff energy scale which we use to calculate other observable manifestations of these 
operators for a number of current and future experimental probes, concentrating on lepton number 
violating phenomena. These include searches for neutrinoless double-beta decay and rare meson, 
lepton, and gauge boson decays. We also consider searches at hadron/lepton collider facilities in 
anticipation of the LHC and the future ILC. We find that some operators are already disfavored 
by current data, while more are ripe to be probed by next-generation experiments. We also find 
. that our current understanding of lepton mixing disfavors a subset of higher dimensional operators. 

' While neutrinoless double-beta decay is the most promising signature of lepton number violation for 

, the majority of operators, a handful is best probed by other means. We argue that a combination 

of constraints from various independent experimental sources will help to pinpoint the "correct" 
model of neutrino mass, or at least aid in narrowing down the set of possibilities. 
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I. INTRODUCTION 



< 

(-H ' The discovery of neutrino masses via their flavor oscillations over long baselines constitutes the flrst solid evidence of 
O . physics beyond the standard model (SM) of particle physics jj. While this is an important first step toward a deeper 
1^ ' understanding of nature, it poses many more questions than it answers. A number of theoretically well-motivated 
models have been proposed and explored to address the origin of the neutrino mass but, strictly speaking, these 
represent only a handful out of an infinite set of possibilities. The question of how well future experiments can probe 
and distinguish different scenarios arises naturally and is quite relevant given the current state of high energy physics. 
The coming years promise detailed explorations of the terascale with the Large Hadron Collider (LHC) and the more 
distant International Linear Collider (ILC) or variants thereof. Expectations are that combined information from 
these two facilities, coupled with high precision, low energy results and cosmological observations will shed light on 
' some of the current mysteries of physics, including that of the neutrino mass. 
: Here, we concentrate on the possibility that the neutrino masses are generated at some high energy scale A where 
T-H U{\) B-L, the only non-anomalous global symmetry of the standard model, is broken. Such a scenario is well-motivated 
, by the observed properties of the light neutrinos including tiny masses, large mixings and the fact that neutrinos are 

■ the only electrically neutral fundamental fermions. More specifically, once U{1)b-l is broken, neutrinos are not 
t — ■ ' protected from getting non-zero Majorana masses after electroweak symmetry breaking. On the other hand, since the 

■ renormalizable minimal standard model* preserves C/(1)b_l, B — L breaking effects will only manifest themselves at 
I low energies through higher dimensional operators. This being the case, one generically expects neutrino masses to be 

suppressed with respect to charged fermion masses by (v/ A)", n > 1, where v is the Higgs boson vacuum expectation 
value. 

By further assuming that all new degrees of freedom are much heavier than the weak scale, we are guaranteed that, 
regardless of the details of the new physics sector, all phenomena below the weak scale are described by irrelevant, 
higher dimensional operators. In this spirit, the observable consequences of all high energy models that lead to 
small Majorana neutrino masses can be catalogued by understanding the consequences of irrelevant operators that 
break B — L hy two units. With this in mind, we will survey all such non-renormalizable effective operators for 
phenomenological signatures at future and current experiments. We restrict ourselves to operators that will lead 
to lepton number violation (LNV), as these will be directly connected to the existence of small Majorana neutrino 
masses. This means that we do not consider operators that conserve L but violate B, and hence also B — L, by 
two units (such operators lead to, for example, neutron-antineutron oscillations), nor do we include operators that 
respect B — L. Most of the time, the latter will not mediate any observable consequences for large enough A, except 
for operators of dimension-six and above that can mediate proton decay. 



* Throughout, we will assume that the weak scale degrees of freedoms are the known standard model fields, plus a minimal Higgs sector. 
Hence, we assume that there are no gauge singlet "right-handed neutrino" fermions or higher SU(2)i^ Higgs boson representations, such 
as Higgs boson triplets. 
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To begin, we systematically name and classify all relevant LNV operators. Fortunately, this has already been 
done^ in up to and including operators of mass dimension eleven. ^ For each operator we then calculate/estimate 
the analytic form of the radiatively generated neutrino mass matrix. Upon setting this expression equal to the 
experimentally measured neutrino masses, we extract the energy scale A associated to the new LNV physics. Armed 
with these scales, we proceed to calculate each operator's phenomenological signatures at a variety of experimental 
settings. Additionally, having explicitly calculated the operator-induced neutrino mass matrices, we may also verify, 
under some generic assumptions, whether one can account for the observable lepton mixing pattern. After such a 
general survey, one is adequately equipped to take a step back and select phenomenologically/theoretically interesting 
operators for further detailed study by "expanding" effective vertices to reveal particular ultraviolet completions. In 
this way, one can use the results presented here as a means of systematically generating renormalizable models with 
well-defined experimental predictions. 

This paper is organized as follows. Sec.|TT]is devoted to an introduction to the effective operators and methods. In 
Sec. Ill Al we derive and comment on the scales A of new physics that are used throughout the remainder of the text. 
In Sec. IIIH we survey various experimental probes of LNV for each operator, and address if and when our analysis 
breaks down due to added model structure or additional assumptions. Specifically, we study both current constraints 
and future prospects for neutrinoless double-beta decay experiments in Sec. IIII Al followed, in Sec. IIII Bl by a similar 
analysis of other rare decay modes, including those of various mesons and W/Z gauge bosons. In Sec. IIII CI we 
present collider signatures of LNV as they apply to future linear collider facilities running in the e~e~ collision mode, 
and describe extensions of our analysis to include associated 77 collisions. We also comment on searches for LNV in 
future hadron machines. Sec. IIVI describes current constraints from neutrino oscillation phenomenology due to the 
general structure of the derived neutrino mass matrices. In Sec. |Vl we highlight a number of "interesting" operators, 
defined by low cutoff scales and prominence of experimental signatures, which are still allowed by current constraints 
on LNV. We undertake a slightly more detailed discussion of their characteristics and signatures and present some 
sample ultraviolet completions. We conclude in Sec. IVII with a summary of our assumptions and results, augmented 
by commentary on future prospects for LNV searches. Our results are tabulated by operator name in Table |T] for easy 
reference. 

We hope that this analysis will prove useful to various audiences on a number of distinct levels. In the most 
superficial sense, the casual reader should note the general features of LNV as well as the diversity of model variations. 
Such information is best expressed in terms of the operator distribution histograms scattered throughout the text. 
These are color-coded by operator dimension or cutoff scale, and typically contain additional information, including 
current experimental prospects. On the more technical side, those interested in specific neutrino mass generating 
models will find detailed, operator specific, information that may be utilized as crude model predictions. Additionally, 
as already alluded to, one may even "hand-pick" operators for model development based on specific phenomenological 
criteria. Finally, we urge experimentalists to search for new physics in all accessible channels. It is our ultimate goal 
to provide motivation for experimental considerations of non-standard LNV effects, beyond neutrinoless double-beta 
decay. 

II. THE LEPTON NUMBER VIOLATING SCALE 

Here we analyze SU{3)c x SU{2)l x U{1)y invariant AL = 2 non-renormalizable effective operators of mass 
dimension up to and including eleven. They are composed of only the SM field content as all other, presumably 
heavy, degrees of freedom are integrated out. As already emphasized, we do not allow for the existence of SM singlet 
states (right-handed neutrinos) or any other "enablers" of renormalizable neutrino masses, such as Higgs SU{2)l 
triplet states. We therefore assume that all lepton number violation originates from new ultraviolet physics and that 
neutrino masses are generated at some order in perturbation theory. 

A d-dimensional operator is suppressed by d — 4 powers of a mass scale A that characterizes the new physics. 



t The authors of 0| discuss all possible effective operators of dimensions up to and including eleven, but only explicitly list those deemed 
unique in the sense that they cannot be written as the product of any previous operator with a Standard Model interaction. We append 
their list and naming scheme to include these into our analysis. 

■t We will argue later that irrelevant operators with mass dimension thirteen and higher, if related to neutrino masses, will require new 
physics below the electroweak scale so that we would have already observed new physics if neutrino masses were generated in this way. 
Furthermore, from a model building perspective, it is difficult to develop models that predominantly yield effective operators of very 
high mass dimension. The probability that such scenarios are both theoretically well-motivated and evade all observations appears to 
be slim. 
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in addition to a dimensionless coupling constant A: 

i 

where we sum over all possible flavor combinations that make up the same "operator-type," as defined below. For each 
operator, A/ A is approximately the maximum energy scale below which the new perturbative ultraviolet physics is 
guaranteed to reside, and A is used as a hard momentum cutoff in the effective field theory. Among all d-dimensional 
operators, we define A so that the largest dimensionless coupling A is equal to unity. Unless otherwise noted, we will 
assume that all other A are of order one. 

In the first two columns of Table HI we exhaustively enumerate all possible lepton number violating operators of 
mass dimension less than or equal to eleven. All together, this amounts to 129 different types of operators, most of 
which, 101 to be exact, are of dimension eleven and consist of six fermion and two Higgs fields. Remaining are 21, 6 
and 1 operator of dimension nine, seven, and five, respectively. The dimension-nine operators can be of two different 
kinds, as defined by their respective field content. They either contain four fermion and three Higgs fields or simply six 
fermion fields with no Higgs field content. For consistency, we use the notation of reference [2], where such a listing was 
first introduced. Our operator naming scheme is also derived from the same list, which we trivially extend to include 
21 elements only mentioned in that analysis. These are the dimension-nine and dimension-eleven LNV operators that 
can be constructed from the "product" of the previously listed dimension-five and dimension-seven operators with the 
SM Yukawa interactions. We individually identify those operators with the same field content but different SU (2)^ 
gauge structure with an additional roman character subscript added onto the original designation from 0. This is 
done in order to render our discussion of the various operators clearer, since specific gauge structures can play an 
important role in the derived energy scale and predictions of a given operator. Note that we neglect effective operators 
that contain SM gauge fields, since, as argued in 0, these are not typically generated by renormalizable models of 
new physics. 

Our notation is as follows. 

are the left-handed lepton and quark SU{2)l doublets, respectively, e"^, u'^ and d° are the charge-conjugate of the 
SU(2)l singlet right-handed charged lepton and quark fermion operators, respectivelty. Conjugate fields are denoted 
with the usual "bar" notation {L, Q, e'^). For simplicity, we are omitting flavor indices, but it is understood that each 
matter fermion field comes in three flavors. All matter flelds defined above are to be understood as flavor eigenstates: 
all SM gauge interactions, including those of the VF-boson, are diagonal. Without loss of generality, we will also define 
the L and e"^ fields so that the charged-lepton Yukawa interactions are flavor-diagonal. 
We take the SU{2)l doublet Higgs scalar to be 

= ( ^0 ) , (11.3) 

and assume that, after electroweak symmetry breaking, its neutral component acquires a vacuum expectation value 
(vev) of magnitude v « 0.174 TeV,* thus spontaneously breaking the electroweak gauge symmetry SU(2)l x U{1)y 
C/(l)cm- In Table m the components of the SU{2)l doublets are explicitly listed and labeled with i,j,k,... = 1,2. 
In order to form gauge singlets, operators are contracted either by the antisymmetric tensor e^ , deflned such that 
£12 — 1, or by trivial contractions with a conjugate doublet field. Different gauge contractions are partially responsible 
for the wide variety of operator structures encountered in this study. 

In order to avoid unnecessarily messy expressions, several features are missing from the operators as listed in Table 
m To begin, SU{3)c color indices are suppressed in these expressions. Color contractions are only implied here 
because SU{3)c is an unbroken symmetry of the SM and hence there is no sense in distinguishing the various quark 
field components. We assume that the parent ultraviolet completion to each operator treats the color gauge symmetry 
properly by introducing appropriately chosen heavy colored particles to render the theory gauge invariant. Slightly 
more serious is the omission of flavor indices to label the fermion generations. For most of this analysis, we assume 
that all new physics effects are generation universal and thus, flavor independent. This is not guaranteed to be the 



* Our numerical value for v is distinct from many treatments of the SM where v is taken to be 0.246 TeV. These are equivalent up to a 
factor of \/2 and are both valid provided a consistent treatment of the interaction Lagrangian. 
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case, as is painfully obvious within the SM. One will also note that, depending on the SU{2) structure of the effective 
operator, different flavor-dependent coefficients will be strictly related. For example, including flavor dependent 
couplings AJ,^, Oi should read \]^ijL]^L-'i^H'^ H'' CikCju where A^^ = A^^^ (symmetric) for all a, /? = e, /x, r. On the other 

hand, O^a should read (for fixed Q and flavors) X^^Ll^L-'^Q'^d'^H^eijeki, where X^^^ = — A^'J^ (antisymmetric) for all 
a,p = e,fi,T. Large differences among the various flavor structures of each operator may very well exist. Flavor is 
an important facet of LNV phenomenology, and is addressed where relevant within the text. 

The final feature missing from our notation is explicit Lorentz structure. Each operator must, of course, form a 
Lorentz scalar, but there are numerous field configurations that can bring this about. The Higgs field is a scalar, and 
as such, transforms trivially under the Lorentz group and is thus of no relevance to this discussion. The fermions, 
however, transform non-trivially and their contractions must be accounted for in each operator. Simple combinatorics 
dictate that there are at most 45 such possibilities for the six-fermion operators that comprise the bulk of our sample, 
3 in the four-fermion case and only 1 for the lone dimension five operator. Additionally, each contraction can be 
made in a variety of ways, corresponding to the bilinear Dirac operators 1, ■y^ and cr'"' = ^[l'^,!''] of the scalar, 
vector and tensor types, respectively. Since we are dealing with chiral fields, the addition of the 75 matrix to form 
the pseudoscalar and axial-vector bilincars is redundant. While this helps reduce the number of possibilities, the task 
of listing, categorizing, and analyzing all possible Lorentz structures for each operator is still quite overwhelming and 
is not undertaken in this general survey. Fortunately, different Lorentz structures for the same operator-type lead to 
the same predictions up to order one effects. This is especially true for the "interesting" operators characterized by 
TeV A scales. We shall quantify this statement and mention specific structures when relevant. That being said, the 
Lorentz structure of an effective operator can suggest a lot of information about its parent renormalizable model. For 
example, it can suggest the spin of the heavy intermediate states and the forms of various vertices. 

Armed with these operators, we can calculate the amplitude of any AL = 2 LNV process. It is important to 
emphasize; that when addressing the phenomenological consequences of any particular operator O, we assume that it 
characterizes the dominant tree-level effect of the new heavy physics, and that all other effects - also characterized by 
other LNV effective operators of lower mass dimension - occur at higher orders in perturbation theory. Our approach 
is purely diagrammatic, in that we begin with an operator-defined vertex and then proceed to close loops and add 
SM interactions as needed to yield the correct external state particles. In this sense, special care must be taken to 
respect the chiral structure as defined by each operator. In order to reach the intended external states, to couple to 
particular gauge bosons, or to close fermion loops, one must often induce a helicity flip with a SM mass insertion. We 
express these inserted fermion masses in terms of the respective Yukawa couplings, yf {f = £,u, d) and the Higgs vev, 
V. The Higgs field can be incorporated into this procedure in a number of ways. We treat the two charged and single 
neutral Nambu-Goldstone Higgs bosons, and respectively, within the Feynman-'t Hooft gauge as propagating 
degrees of freedom with electroweak scale masses. The physical neutral Higgs, ho, can be either chosen to propagate 
as a virtual intermediate state, or couple to the vacuum with amplitude v. 

In order to avoid the task of explicitly evaluating a huge number of multiloop Feynnman diagrams, we succumb 
to approximate LNV amplitudes based on reasonable assumptions and well-motivated rules. Our methodology is 
motivated by exact computations with one-loop, dimension 7 operators where the work is analytically tractable, as 
well as on general theoretical grounds. For select operators, we have also checked our assumptions against predictions 
from ultraviolet complete models with success. In order to perform a particular calculation, we draw the appropriate 
diagram(s), taking care that no momentum loop integral vanishes by symmetry reasons. This step is potentially 
quite involved, as multiple diagrams can give sizable amplitude contributions depending on the characteristic energy 
transfer in the system, not to mention the cumbersome Dirac algebra within the respective loops. Given the high, 
often super- TeV, mass scale associated with our calculations, it is often convenient to work in the gauge field basis 
where each boson state is associated with a single SM group generator, as is natural before electroweak symmetry 
breaking. In a similar sense, all fermions, iuciluding those of the third generation, are taken to be massless to zeroth 
order. All masses are included perturbativcly where needed via mass insertions. At first guess, it would seem that 
our results are only valid in the rather subjective limit A ^ By direct comparison with other more complete 
approximations, however, we find that our predictions are very reasonable at all scales above 0.5 TeV. Keeping all of 
this in mind, we apply the following "rules" to obtain approximate amplitude expressions. 

1. Trivial numerical factors: A number of numerical factors can be read off trivially from the Feynman diagrams. 
Specifically, one can extract the presence of the suppression scale A^'^'^^^' directly from the dimension d operator, 
as well as the dimensionless coupling constants A. Generally, A is a generation dependent quantity, but for lack 
of any experimental evidence to the contrary, we take A = 1 universally unless stated otherwise. In the case of 
scenarios already constrained by current data, we will relax this assumption to "save" the operator and comment 
on the phenomenological consequences of the change. 

Furthermore, various factors of the electroweak scale v may be extracted from the operator's Higgs field content, 
in addition to fermion/gauge boson mass terms. In this way, we may also include the various Yukawa and gauge 
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coupling factors Vf [f = i,u,d) and gi, respectively, where i runs over the three SM gauge groups. For 
simplicity, we neglect the gauge subscript i in further analytic expressions. Finally, a color factor of 3 associated 
with each quark loop should also be included in our computations, but can (and will) be neglected for simplicity 
from algebraic expressions where order one factors are irrelevant and only serve to render expressions more 
cumbersome. We note that all coupling constants are subject to renormalization group running. In particular, 
those occurring within a loop should be evaluated at the scale A. We neglect this order one effect since it is 
most important at large A scales where operators tend to have less of a phenomenological impact due to the 
(1/A)" suppression. 

2. Loop factors: In all of our calculations, we assume that each operator defines an effective field theory, charac- 
terized by the scale A. This implies that all momentum integrals are effectively cut off at A, above which new 
states will emerge to regularize the theory. Divergences in such loops tend to cancel the large scale suppressions 
inherent to the bare operators, and thus enhance predicted LNV rates. Specific divergences can be determined 
by simple power counting of momentum factors. Of course, multiple loop integrals are often convoluted to the 
point where substantial simplification is needed to determine the dominant divergent term. Such a complication 
is in part due to the numerator of the Dirac propagators, which include single momentum factors and must 
therefore be present in pairs to contribute effectively to an ultraviolet divergence. The process of adding loops 
to induce A power law divergences should only be pursued to the point where the suppression of the induced 
effective term is no less than A~^. Any further divergent contribution must be treated as a renormalization to 
lower order terms, and hence, can only add small finite corrections to the total amplitude. In any case, those 
diagrams with the smallest scale suppressions are not always the most dominant, as will become clear later when 
we discuss specific results. 

In addition to power-law divergences, each loop is also associated with a numerical suppression factor. This 
arises from the proper normalization of the loop four- momentum integral as a factor of (27r)^^, the characteristic 
phase space "volume" of a quantum state. It allows one to view the integral as a coherent sum over all possible 
intermediate configurations in a consistent way. Partially evaluating these integrals for a number of examples, 
one quickly finds that two powers of tt cancel with the four dimensional Euclidian space solid angle / dfl4. We 
introduce a suppression factor of (167r^)~^ 0.0063 for each diagram loop, which tends to offset enhancements 
from associated divergent factors. A quadratically divergent loop diagram is often proportional to the lowest 
order contribution times (l/167r^)(A/u)^ to the power n (number of loops in the diagram). This contribution 
is larger than the leading order one if A > inv ^ 2 TeV for any number of loops. The situation is often more 
involved, as many loops turn out to be logarithmically divergent or even convergent. The important conclusion is 
that adding loops is not an efficient way to enhance LNV rates at the low scales accessible to future experiments. 
This fact is demonstrated by example in Sec. IIIII 

Finally, as already alluded to, many diagram loops will exhibit logarithmic divergences, as is the standard case 
in renormalizable theories involving fermion and vector fields. This occurrence typically reflects the differences 
between the two characteristic scales inherent to the system, namely A and v, and are of the general form 
^log"(A/w) for some power n. ^ is a small, loop suppressed, dimensionless coupling coefficient. Numerically, 
these logarithms are much softer than their quadratically divergent counterparts seen elsewhere in the diagrams 
and can safely be neglected. 

3. Intermediate states: We treat all virtual intermediate states, outside of loops, as if they carry the characteristic 
momentum of the interaction Q and neglect Dirac structure, unless stated otherwise. In particular, goldstone 
bosons are assigned the propagator {Q^ — Mg)~^ and fermions are assigned (Q — Mf)^^. In the case of an 
intermediate neutrino, this reduces to a simple factor of for all realistic Q values. Hence, for very low 
energy processes {Q ^ 100 MeV), neutrino exchange diagrams tend to dominate LNV rates. 

4. Lorentz structure: For the purposes of our analysis, we assume that all Lorentz contractions between fermions 
are scalar-like. As previously mentioned, the absolute magnitude of most LNV amplitudes is robust under 
this assumption up to order one factors. The only qualitative exception to this occurs in some cases involving 
fermion bilinear terms with a tensor Lorentz structure {ip'cTp^i/^')- This factor, when coupled between two fermions 
contracted in a loop, will yield a vanishing rate due to its antisymmetry inside of a trace, since Tr^a^^^) — 0. 
This can be bypassed by introducing a new momentum vector into the trace, implying the addition of another 
loop. In most cases, this is most efficiently accomplished with a new gauge boson line, which is accompanied 
by a logarithmic divergence. The combination of both factors leads to a marginal amplitude suppression (with 
respect to the same operator where all fermion bilinears are Lorentz scalars) for all energies of interest. 

With these approximations in hand, it is a simple matter to estimate the amplitude associated to any given diagram. 
Still, one must wonder about the uncertainty induced onto the calculations by such varied assumptions. Can results 
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obtained by such methods supply vahd physical predictions? The answer, of course, depends on the question that 
is being asked. Here, we will only be interested in estimating order of magnitude effects, including what value of 
A is required in order to explain the observed neutrino masses and, once A is so constrained, what is the order of 
magnitude of other related observable effects. 

One may wonder whether a more detailed estimate of the effects of each individual operator would lead to more 
reliable results. The answer is negative. It is easy to show that different renormalizable theories that lead to the 
same effective operator at tree-level will mediate different processes at the loop-level with order one different relative 
strengths. Furthermore, the derived cutoff scales inherit the uncertainty from the absolute value of the heaviest 
neutrino mass, which is only loosely bounded between 0.05 eV and 1 eV by the extracted atmospheric mass squared 
difference @ and tritium beta decay kinematic measurements 0, [E| • This is an order of magnitude uncertainty that 
cannot be avoided even if one were to perform a detailed computation within a well-defined ultraviolet complete 
theory. 

In summary, given all approximations and uncertainties, our results are only valid up to ± an order of magnitude. 
In this spirit, one need not explicitly consider order one factors that will necessarily yield negligible corrections by 
these standards. Such a large error tolerance supplies the need for care when interpreting results. In particular, one 
should not place too much emphasis on any one bound or prediction, unless it is very robust, i.e., able to withstand 
variations of at least a factor of ten. Of course, for those operators constrained by several different independent 
sources one can, and should, take more marginal results seriously. 



A. Neutrino Masses and the Scale of New Physics 



Having defined the set of LNV operators, we now extract the scale of new physics from the direct comparison of 
radiatively generated neutrino mass expressions to their observed values. Since there are three light neutrino masses, 
we will use the heaviest of these to set the overall mass scale. Neutrino oscillation data, currently providing the only 
evidence for neutrino masses, constrain the relative magnitudes of the mass eigenstates but not the overall scale jlj. 
Such data only supply a lower bound on the heaviest neutrino mass, derived from the largest observed mass squared 
difference Am^g = |m§ — mf \ w 0.0025 eV^, the atmospheric mass squared difference ^Sjj. At least one neutrino mass 
must be greater than Am^g « 0.05 eV. Neutrino oscillations also teach us that the next-to-heaviest neutrino weighs 
at least i/Am^ « 0.009 eV (the solar mass-squared difference), in such a way that the ratio of the heaviest to the 
next-to-heaviest neutrino masses is guaranteed to be larger than, approximately, 0.2. No lower bounds can be placed 
on the lightest neutrino mass. An upper bound on the heaviest neutrino mass is provided by several non-oscillation 
neutrino probes. Cosmology provides interesting constraints on the sum of light neutrino masses, but these are quite 
dependent on unconfirmed details of the thermal history of the universe and its composition [0]. Most direct are 
kinematic measurements of the tritium beta decay electron endpoint spectrum [t^. Both types of probes provide 
upper bounds near 1 eV, likely to improve in coming years. We choose to perform our calculations assuming the 
mass scale TOj/ ~ 0.05 eV, corresponding to the experimental lower bound. In this way, each extracted operator scale 
A, inversely related to the neutrino mass, represents a loose upper bound. Since most rates for LNV observables are 
proportional to some inverse power of A, this choice implies the added interpretation that, all else remaining equal, 
our results for such rates should conservatively reflect lower limit predictions. 

LNV neutrino masses are nothing more than self-energy diagrams evaluated at vanishing momentum transfer. 
These must couple together the left-handed neutrino state Va with the right-handed anti- neutrino state Vfs, as shown 
schematically in diagram (a) of Fig. [T] Here the flavor indices a and /3 can accommodate any of the three lepton 
flavors (a, (3 = e, ii,t). The derived Majorana masses TOq/j = mpa are generally complex. The large grey circle in this 
diagram represents all possible contributions to the neutrino mass. Specifically, it contains the underlying AL — 2 
operator along with all modifications needed to yield the correct external state structure. This includes such objects 
as loops, additional gauge boson propagators and SM coupling constants. Generally, several diagrams can contribute 
to this mass generation, but special care must be taken that these are not proportional to any positive power of Q, 
the momentum carried by the neutrino legs, as this would not lead to a nonzero rest mass correction. 

Diagrams (6) — (e) of Fig. [T] are examples that serve to illustrate some typical features encountered in our effective 
operator induced self-energy calculations. The underlying LNV operators shown in each diagram contain six fermion 
fields and are therefore of dimension nine or eleven, as are the majority of the analyzed operators. Each of these 
diagrams generates an effective dimension-five interaction 

_ (5)(i^ff)(L^ 
^5-4a/3 : 

where ^^^^ = is a generation dependent coupling constant that is calculable, given the structure of the original 
operator. This is easily verifiable via direct power counting, despite differences in dimension, loop number, field 
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(a) 



FIG. 1: Sample diagrams that radiatively generate Majorana neutrino masses. Diagram (a) is representative of all operators 
that can generate the needed external state neutrinos. This usually proceeds via loop contractions and other couplings, hidden 
within the light gray region. Diagrams (6) — (e) help illustrate the methodology of this analysis. Despite obvious differences, 
all of these generate effective dimension-five interactions of calculable strength. See the text for more details. 



content, and helicity structure. It turns out that most operators, especially those characterized by super-TeV scales, 
possess this property. 

We describe each sample diagram in turn to point out important features. A subset of the subtleties described 
below is encountered when estimating the neutrino masses map for the entire effective operator set. Diagram (6) 
is a simple two-loop radiatively generated mass term proceeding from dimension-nine operators, such as Oiij, = 
U U Q'^ d'^Q'^ d'^eik^ji , containing the fermion structure fhfiif'LfR: where the fields / and /' are contracted into loops 
with mass insertions that supply the needed field chirality flip. Masses arising from such operators are proportional 
to two powers of fermion Yukawa couplings. Strictly speaking, allowed fermions from all three generations traverse 
the closed loops and contribute to the mass. However, assuming universal new physics coupling constants, third 
generation fermions will strongly dominate the induced neutrino mass. In cases such as these, one can freely suppress 
couplings to the lighter two generations without modifying the expected value of A. Since diagrams arising from the 
majority of our operator set contain at least one loop of this kind, this property proves quite useful when attempting 
to avoid low energy nuclear physics constraints, as will be discussed in more detail later. 

Diagram (d) involves an operator of dimension eleven, such as O22 = L'-U L''e'^Lke'^H''H™eiiejm, but has a similar 
structure to Diagram (6) since both neutral Higgs fields Hq couple to the vacuum, yielding a v'^ factor. In this case, 
the parent underlying operators contain the fermion structure fL{R)fi^(ii)fL(^if')f'£(^ijp or simple variants thereof. From 
this, it is clear that such operators will create and annihilate the same field and one can close the fermion loops 
without mass insertions. A little thought reveals that such loops, if left on their own, will vanish by symmetry, 
since J d'^k[k^ /g{k'^)] — for all functions g{k'^). Hence, non-zero neutrino masses only appear at a higher order in 
perturbation theory (i.e., we need to add another loop). To maintain the chiral structure of the diagram a gauge 
boson line insertion is always the most effective. The specific gauge field required in this step depends critically on the 
quantum numbers of the fermions / and /' contained in the operator itself and must be determined on a case-by-case 
basis. The absence of Yukawa dependence renders the estimated value of the cutoff scale A insensitive to the values 
of the dimensionless operator couplings (A), given the way A is defined. Notice that this three- loop diagram, like 
Diagram (6), predicts an anarchic neutrino Majorana mass matrix, currently allowed by the neutrino oscillation data 
0. That is, up to order one corrections, all entries are of the same magnitude, nii, « 0.05 eV. This is in contrast to 
the remaining sample diagrams ((c) and (e)), which both suggest flavor-structured mass matrices. 

Dimension-nine operators, such as Oig = L^Q^ d'^d'^e'^u'^eij , yielding diagram (c) have the peculiar property that, 
upon expanding out the various SU{2)l contractions in terms of component fields, no Vai^p content is present to form 
the external legs of a mass diagram. Here the LNV is introduced via the fermion structure Va&R.p, which annihilates a 
left-handed neutrino and creates a left-handed positron. Hence, to tie in the needed antineutrino line, one must both 
flip the charge lepton helicity and carry away the excess charge with some bosonic state. Of course, such a charged 
boson is guaranteed by charge conservation to be needed elsewhere in the system to close some ///' loop. In this 
particular example, the process is illustrated by the exchange of a charged Higgs goldstone boson H~ . The crucial 
point is that this mass is necessarily proportional to a charged lepton Yukawa coupling y^^ , of the same flavor as 
the external neutrino, since we are working in the weak eigenbasis where the gauge couplings and the charged-lepton 
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Yukawa couplings are flavor diagonal. By symmetry, the contribution to the mai3 entry of the neutrino mass matrix is 
proportional to yi^ + yi^ , which reduces to the largest coupling y^^ in the realistic case of hierarchial charged lepton 
Yukawa couplings. 

Finally, Diagram (e) yields a five-loop suppressed neutrino self-energy originating from a dimension-eleven LNV 
operator such as Osg = e'^e^Q'^d'^Q^ dl^H^ H^eikCji. This represents the most complicated structure considered in this 
analysis. As in Diagram (c), no explicit VaVjs structure is available in the underlying operator, but in this case all 
of the LNV arises from e^^-^e^^-type interactions. Curiously, this interaction already flips helicity as it annihilates a 
left-handed positron and creates a right-handed electron. Unfortunately, being an SU{2)l singlet, Cfl only couples to 
the neutrino via a charged Higgs induced Yukawa interaction; therefore, this amplitude must be proportional to the 
product yi^yii^ to yield a legitimate neutrino mass contribution. One might imagine that the Higgs fields contained 
in the LNV operator could be used to produce the needed neutrino legs, but this is not possible since the resulting 
loop would have a structure of the form J c?'*fc-^^^tS_ which vanishes in the "rest mass" limit. It is clear that 

both Higgs fields must again couple to the vacuum and the needed flip must come from the other fermion loops. The 
resulting loop integrals can be separated into two convoluted pieces corresponding to both loop/leg pairs. A little 
thought reveals that each loop set contains three fermion lines whose associated integral is again proportional to the 
momentum of the external neutrino, and thus is not a valid mass correction. To fix this last problem without further 
complicating the chiral structure, one can add a gauge boson exchange between the fermion loops, as was also done 
in Diagram (d). 

Despite the dominance of the generated dimension-five interactions described by Eq. (jll.4p for the majority of the 
studied LNV operators, we find that this need not be the case for all of them. For some operators, the dimension-five 
neutrino-mass effective operator Eq. pi.4[) occurs at higher order in perturbation theory than the dimension-seven 
neutrino-mass effective operator (schematically, {LH)^H^). For these, the neutrino mass matrix is generated after 
electroweak symmetry breaking from 

{LH){LH) (7) {LH)(LH)HH 
^" = 16^ A + A5 ' ^"-^^ 

(7) 

where are new calculable coefficients. This type of structure is present in the following operators 

O7, 021^,, , O22, O23, 025,026, , O27, , 029^030, , 031,044. , O57. (H.6) 

In general, they are associated with dimension-eleven operators^ whose mass diagrams are found trivially by connecting 
the external fermion loops and coupling the neutral Higgs fields to the vacuum. This adds two factors of the electroweak 
scale to the mass expressions. Dimensional analysis dictates that the fermion loops must conspire to yield an addition 
factor of , usually from mass insertions utilized to flip helicities. For the dimension-seven operators in Eq. pi.5[) . 
the resulting neutrino mass expression is proportional to w^/A^. If we assume, as is usually the case, that most of the 
dimensionless factors of Eq. (jll.sp are common to both and we find rriap oc l/167r^ -f j . In such cases, 
the dimension-seven contribution is only relevant for operator cutoff scales A < 47rw « 2 TeV. Such low scales are 
seldom reached considering that these operators are efficient at mass generation at low orders and consequently do 
not possess the necessary suppression factors. Still, for completeness, we include these terms when relevant. 



TABLE \: Dimension-five through dimension-eleven LNV operators analyzed in this survey. The first two columns display the 
operator name and field structure, respectively. Column three presents the induced neutrino mass expressions, followed by 
the inferred scale of new physics, A^. Column five lists favorable modes of experimental exploration. Column six describes an 
operator's current status according to the key U (Unconstrained), C (Constrained) and D (Disfavored). See text for details. 



o 


Operator 




A. (TeV) 


Best Probed 


Disfavored 


1 

2 

3a 
36 


UL^Q''d''H'e,k€ji 


Vli u 
16^ A 

Vdg D 

(leir^)^ A 
Vd v'^ 
TGtv^ a 


6 X 10" 
4 X 10^ 
2 X 10^ 
1 X 10** 


I3I30U 


U 
U 

u 
u 



t This also occurs with operator Oy, which is of dimension nine. This operator is the exception, in that it explicitly contains three Higgs 
bosons which naturally aids in building the needed v'^ factors in a way similar to that discussed in the text. 
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TABLE I: Dimension-five through dimension-eleven LNV operators analyzed in this survey. The first two columns display the 
operator name and field structure, respectively. Column three presents the induced neutrino mass expressions, followed by 

the inferred scale of new physics, A^. Column five lists favorable modes of experimental exploration. Column six describes an 
operator's current status according to the key U (Unconstrained), C (Constrained) and D (Disfavored). See text for details. 



o 


Operator 


mal3 


K (TeV) 


Best Probed 


Disfavored 






16^ A 




4 X 10'^ 




U 






yu9 v'' 

(16ir^)^ A 




6 X lO'' 




u 


5 




Vd 

(ISir^)^ A 




6 X 10^ 


/3/?0z/ 


u 


6 




(16ir^)^ A 




2 X lO'' 


/3/30z/ 


u 


7 




y^fi (16ir^)^ A I 16^ + ] 


4 X 10^ 


mix 


c 


8 




VdVu 




6 X lO''' 


mix 


c 




(levr^)^ A 




9 


UUL''e''L'-e''tijtki 


2 Q 
(16ir^)^ A 




3 X 10^ 


/3/30z/ 


u 


10 


rUL''e''Q'd''eijeki 


VlVd 
(167v^y-^ A 




6 X 10^ 


/3/?0z/ 


u 


11a 


rUQ''d''Q^d''€ij€ki 


VdS u 
(16ir^)a A 




30 


/3/?0z/ 


u 


Hi, 


rUQ''d''Q'd''eikeji 


^,2 2 
(16ir2)2 A 




2 X lO"* 


/3/30i/ 


u 


12a 




u2 .,2 
(16ir^J^ A 




2 X 10^ 


/3/30i/ 


u 


126 




1.2 „2 2 
(lfi7|2)a A 




4 X lO"* 


/3/30j^ 


u 


13 




(167r2)2 A 




2 X 10^ 


/3/30z/ 


u 


14a 


UUQi^u'=Q''d''€ij 


(167r2)a A 




1 X 10^ 


/3/30z/ 


u 


146 


DUQiU'^Q^d'^eji 


VdV^ 
(167r2)2 A 




6 X 10^ 


/3/30z/ 


u 


15 


U'UL'^d''TiU<'ejk 


VdVug'' 

(16ir2)3 A 




1 X 10^ 


/3/30i/ 


u 


16 




(16ir2)4 A 




2 


/3/30i/, LHC 


u 


17 


rL^d''d''d'=u'=eij 


4 9 

VdVug D 

(16ir2)4 A 




2 


I3I3QV, LHC 


u 


18 


L V d u u'^u'^eij 


VdVug'^ 

(167r2)4 A 




2 


/3/?0z/, LHC 


u 


19 


T 2 7C TC ~f ~f 

L Q-i d d e'^u'^eij 






1 


ppOz/, HEmv, LHC, mix 


c 


f^/S (167r2)3 A 

Vdvl v'^ 




20 






40 


/3f30u, mix 


c 


(I6,r2-|a A 




21a 


r i T 1 T k C C TTTtl TTTl 

L U L e Q u H H eijekm^n 


(16ir2)2 A I 16^ + 




2 X lO"* 


ppOl' 


u 


2l6 


L''e'=Q''-a'H"'H"eiiejmek„ 


y^.yu v'^ ( 1 1 
(16ir2)2 A 1^16^ ^ 


\ 


2 X 10'^ 


/3/30j/ 


u 


22 


UUL''e''Tke''H'-H"'euejm 


j,2 

(16ir2)a A 




4 X lO"* 




u 


23 


L'L'L''e''Qkd''H^H'^tritj^ 


ViVd V 1 1 1 
(16ir2)2 A I 16^ 


\ 


40 




u 


24a 


UUQ''d''Q'-d''H^Hiejkeim 


.,2 2 
(16ir2)3 A 




1 X 10^ 


/3/30z/ 


u 


246 


UUQ''d''Q'-d''H^Hiejmeki 


^.2 2 
fd w 
(167r2)3 A 




1 X 10^ 


/3/30z/ 


u 


25 


UUQ''d''Q'u''H"'H"€imejn€kl 


UrtVll iP" { 1 1 

if a, a V 1 A 1 

(i67r2)^ A ^^16;;^ 


A^J 


4 X 10^ 


/3/?0z/ 


u 


26a 




yiyd 

(16ir2)a A 




40 


/3I301' 


u 


266 


VUQ''d''Lke-H'H"'eaejm 


yevd / 1 1 

(16ir2)2 A I 16^ ^ 


.2\ 


40 




u 


27a 


VL'Q''d''Q,d-H'H"'eji€km 


9^ «^ 
(16ir2)a A 




4 X lO"* 


f3f30v 


u 


276 


L'L'Q''d''Q^ d'^H^ H'^eu 


2 2 
9 I' 
(16ir2)a A 




4 X 10* 


/3/30z/ 


u 


28a 


UUQ''d''Q^u''H^Hieki 


(167r^)3 A 




4 X 10^ 


/3/30z/ 


u 


286 


UUQ''d''Q^u''H^'Hieji 


ydy-u 

(16ir2)a A 




4 X 10^ 




u 


28c 


VL^Q''d''QiU''H'^Hiejk 


!/dJ/t» 
(16ir2)a A 




4 X 10'^ 




u 


29a 




I/,^ ( 1 1 
(167r^)2 A I 16^ 




2 X 10^ 




u 


296 


L*UQ''u''QiU''H'H"'c,kCj,n 


2 2 

9 

(167r2)3 A 




4 X 10"* 




u 


30a 


L'ULii'=QkU'=H''H''eji 


(167r2)a A 




2 X 10^ 


/3/30z/ 


u 


306 


VL^Lme-Q„u-H''H^eikejie"'" 


yeyu ( i , 

(16ir2)2 A I 16^ ^ 


*1 


2 X lO'"* 




u 


31a 


rUQid'^QkU''H''H^€ji 


!/d!/u f 1 1 
(167r2)2 A \^T6^ 


w ) 


4 X 10^ 


/3/30z/ 


u 
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TABLE I: Dimension-five through dimension-eleven LNV operators analyzed in this survey. The first two columns display the 
operator name and field structure, respectively. Column three presents the induced neutrino mass expressions, followed by 
the inferred scale of new physics, A,^. Column five lists favorable modes of experimental exploration. Column six describes an 
operator's current status according to the key U (Unconstrained), C (Constrained) and D (Disfavored). See text for details. 



o 


Operator 




A„ (TeV) 


Best Probed 


Disfavored 






VdVu I 1 1 u \ 
(16ir^)'-^ A \T6^ ~^ W J 


4 X 10^ 






U 


32a 




,,2 2 
(16ir2)a A 


2 X 10^ 






U 


326 




,,2 2 
(16ir2)a A 


2 X 10^ 






u 


33 




2 2 
(16ir2)a A 


4 X 10" 






u 


34 




y^fi (167r2J<l A 


< 0.5 


f3/30iy, mix, ILC, LHC 


c 


35 


Ue''QjU''WH''e,k 


2 2 
y'fl (167r2)4 A 


2 




mix, LHC 


c 


36 


e'=e'=Q'd''Q^d''H''H't,keji 


yi^yifi (i6ta)B A 


< 0.5 




mix, HElnv, ILC, LHC 


D 


37 


e'=e''Q'd''Q^u''H^H''e,k 


2 2 

VrlVllQ V 


< 0.5 




mi^ HRInv TT P T HP 


D 


38 




Vla^Vlf} (16^^)6 A 


< 0.5 


ppw, 


miY MTTlm^ TT P T PTP 
iiiix, njzjiiiv, 1 j ri v_ 


D 


39a 




2 1 
9 " 
(16ir2)a A 


8 X lO"* 






U 


396 


U U L} LmLnH"^ eijEki 


2 2 
9 " 
(16ir2)a A 


4 X lO"* 






u 


39c 


U U L} LiLmH"^ ejk^in 


2 2 
9 " 
(16ir2)a A 


4 X lO"* 






u 


39d 


L'^ LpLqH"^ H^'eijej^^ein^^^ 


2 2 
9 ^ 

^ 2 ^ 2 


4 X lO"* 






u 


40a 


U U L'' LiQ e^rnUn 


9 ^ 
(16ir2)a A 


4 X lO"* 






u 


40b 




2 2 
9 V 

(16ir2)a A 


4 X 10'' 




RRf\ij 


u 


40c 




2 2 
9 I) 
(16ir2)a A 


4 X lO"* 




ppKjiy 


u 


40d 




2 2 
9 I) 
(16ir2)a A 


4 X 10* 




RRV\i/ 
ppuiy 


u 


40e 




(16ir2)a A 


4 X 10* 






u 


40/ 




^ 2 ' 2 


4 X lO* 




RRO1J 
ppuiy 


u 


4O9 




9 " 
(16ir2)a A 


4 X lO* 






u 


40h 




2 2 
9 " 
(16ir2)a A 


4 X lO* 




RROi^ 
ppuiy 


u 


40i 




2 2 
9 " 
(16ir2)a A 


4 X 10* 




f3f30u 


u 


40i 


rU L''Q'L^Q„H^H'>tiptiqejke"'" 


2 2 

9 

(16ir2)a A 


4 X 10* 




RROij 
ppKjiy 


u 


41a 


rUL''d''L,d-H'H"'ejiekm 


2 2 

9 

(16ir2)a A 


4 X 10* 




f3(30u 


u 


4l6 


U UL''d''Li d'^H'-H"' uj ekrn 


2 9 

9 

2 ' 2 


4 X 10* 






u 


42a 




9 I) 
2 ' 2 


4 X 10* 






u 


426 


L' U L'^u'^Liu'^H'^ H'^eijekrn 


9 I) 
(16ir2)a A 


4 X 10* 






u 


43a 


L'UL''d''Liu''H'Hiejk 


2 9 

VdVug u 

(167r2)<l A 


6 




f3(30u, LHC 


u 


436 


UUL'^d'^TjU'^H'^'Hieki 


ydyug^ 

(167r2)<l A 


6 




f3(30u, LHC 


u 


43c 


UUL''d''Liu''H"'Tir,e,iekme''^ 


Vdy^g^ 

(167r2)4 A 


6 




PPOu, LHC 


u 


44a 


UUQ''e''Q^e''H'H'^ejiekm 


9' «^ 
(16ir2)a A 


4 X 10* 






u 


446 




9' 

(16ir2)a A 


4 X 10* 






u 


44c 




4 2 

9 

(167r2)4 A 


60 






u 


44d 




2 2 
9 ^ 
(16ir2)3 A 


4 X 10* 






u 


45 


L'Ue''d''e'=d''H''H'e,k£ji 


2 2 
9 I) 
(16ir2)a A 


4 X 10* 






u 


46 


L'L'e''u''e''u''H''H'e,ke,i 


2 2 
9 I) 
(16ir2)a A 


4 X 10* 






u 


47a 




9' «^ 
(16ir2)a A 


4 X 10* 






u 



t This operator is modified slightly from its original form as given in reference 3 where it appeared as 039(0) = 
V- LiLjH™-H"ejm<^],l. We corrected this error. 
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TABLE I: Dimension-five through dimension-eleven LNV operators analyzed in this survey. The first two columns display the 
operator name and field structure, respectively. Column three presents the induced neutrino mass expressions, followed by 

the inferred scale of new physics, A^. Column five lists favorable modes of experimental exploration. Column six describes an 
operator's current status according to the key U (Unconstrained), C (Constrained) and D (Disfavored). See text for details. 



o 


Operator 




K (TeV) 


Best Probed 


Disfav 


476 






4 X lO"* 




U 


47, 




2 2 


4 X 10"* 


/3/?0z/ 


u 


47d 




2 2 


4 X 10" 


/3/30z/ 


u 


47e 




2 2 
9 V 

(16772)-^ A 


4 X 10* 


/3/30z/ 


u 


47/ 




4 9 

n67r2)4 ^ 


60 


,3/30z/ 


u 


47« 




2 2 
9 " 
ng7^2^c( A 


4 X lO"* 


/J/JOf 


u 


47h 




4 2 
9 w 

(leir^)* A 


60 


/3/30z/ 


u 


47> 




9^ 

(16-7r2)a A 


4 X lO'' 


/3/?0z/ 


u 


47, 




9^ 

(167r2)3 A 
2 ' 2 


4 X 10* 


/3/30z/ 


u 


48 


U'Ud''(r-d''d''H''H''e,keji 


9 V 
(16ir2)3 A 


4 X 10* 


/3/30z/ 


u 


49 


UUd''u''d''u''H''H''e,keji 


9^ 

(16ir2)3 A 


4 X 10* 


/3/30i/ 


u 


50 


L'Ud''d''d''u''H''~H,t,k 


(16ir2)4 A 


6 


/3/30i/ LHC 


u 


51 




2 2 

9 

(167r2)a A 


4 X 10* 


/3/30z/ 


u 


52 




(16ir2)4 A 


6 


/3/30Z/, LHC 


u 


53 


UUd'=d'=u<=u<=HiHj 


(167r2)5 A 


< 0.5 


/3/30i/, HElnv, ILC, LHC 


D 


54a 


rQ'Q''d''Q^e-H'H"'eji€km 


(167r2)4 A 


< 0.5 


PPOv, mix, HElnv, ILC, LHC 


D 


546 


rQ'Q''d''Q.e-H'H"'euekm 


Vd9^ 
y«3 (167r2)4 A 


< 0.5 


PHOu, mix, HElnv, ILC, LHC 


D 


54c 


rQ^Q''d''Qie<^H^H"'eimejk 


2 2 
Vda J) 
'^'^^ (167r2)<l A 


< 0.5 


PPQv, mix, ILC, LHC 


D 


54d 

55a 


UQiQ''d''Qie''H'-H^tijekm. 


»d9^ 
y*|3 (167r2)4 A 

(167r2)4 A 


< 0.5 
2 


/3/?0z/,mix, HElnv, ILC, LHC 
PPQv, mix, LHC 


D 

C 


556 




2 2 
,,. !/u9 w 
»*/3 (1677^)* A 


2 


/?/30z/, mix, LHC 


C 


55c 




(ISTT^y A 


2 


l3l3Qu, mix, LHC 


C 


56 
57 


UQWd''e''d''H*'H'-tikeji 
L'd''QjU^e''d^H^H''€ik 


y'^H (167r2)4 A 
(167r2^4 A 


< 0.5 
2 


/3/30i^, mix, ILC, LHC 
/3/?0z/, mix, LHC 


c 
c 


58 




9u9 

»''/3 (167r2)4 A 


2 


mix, LHC 


c 


59 


V d''d''e''u''H''Tli ejk 


(167^2)4 A 

VdVu 
{167r2)4 A 


< 0.5 


PPOv, mix, HElnv, ILC, LHC 


D 


60 


rd^QjU^e^u^H^Hi 


< 0.5 


PPQv, mix, HElnv, ILC, LHC 


D 


61 




2 / , 2 \ 
167r2 A I 167r2 A^ J 


2 X lO '' 


I3I3QV 


u 


62 


L'UL''e''H'L''e''Hrtijtki 


^,2 2 / , 2 \ 
V 1 I 1^1 
(167r2)2 A \T6^ ^ 7^ ) 


20 




u 


63a 


L^UQ^d^H^Ue^Hreijeki 


y^Vd v'^ 

(167r2)a A 


40 


/3/?0z/ 


u 


636 


T i T i /^fc TC TtI T V C TT 

L U Q d H L' e Hr€ikeji 


ytvd ( 1 1 \ 

(167r2)2 A IT6^+ A^i 


40 




u 


64a 


r i T i ~f Tjk T T C TT 

L L-' Q^u'^H L e Hrejk 


VtVu ( 1 1 I 
(167r2)2 A "t A^ y 


2 X 10^ 


ppQv 


u 


046 




veVu v'^ 

(167r2)a A 


Owl n3 
Z X iU 




TT 

u 


65 


Ve''u''d''WUe''lirtr] 


VdVug'^ 

(167r2)4 A 


6 


Pl30v, LHC 


u 


66 


UUH''H'tikQ^d''lir^,i 


Vd ( 1 1 ^ 
167r2 A 167r2 A^ ^ 


6 X 10^ 


I3l30u 


u 


67 


L'UL''e''H''Q''d''Hre^jeki 


yiVd ( 1 \ 
(167r2)2 A ^ ) 


40 


/3/30v 


u 


68a 


L'L'Q''d''H'Q'-d''Hreijeki 


(167r2)a A ^ 1? ) 


1 


Pl30u, LHC 


u 


686 


U UQ^d^H^Q'-d'^H reiki ji 


""l ( 1 1 
(leir^)^ A \1E^^J?] 


1 X 10^ 


/3/30z/ 


u 


69a 


DUQiU'^H'^Q'-d'' Hrejk 


VdVu ( 1 1 1 
(167r2)2 A \T6^ ^ 7^ ] 


4 X 10^ 


/3/?0z/ 


u 


696 


UUQ^u'=H^Q''d''Hreij 


VdVua^ v'^ ( 1 1 f ^ 1 

(167r^)a A ^^16;;^ A^-y 


7 


/3/30Z/, LHC 


u 
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TABLE I: Dimension-five through dimension-eleven LNV operators analyzed in this survey. The first two columns display the 
operator name and field structure, respectively. Column three presents the induced neutrino mass expressions, followed by 
the inferred scale of new physics, A,^. Column five lists favorable modes of experimental exploration. Column six describes an 
operator's current status according to the key U (Unconstrained), C (Constrained) and D (Disfavored). See text for details. 



o 


Operator 




K (TeV) 


Best Probed 


Disfavored 


70 
71 
72 

73a 

736 

74a 

746 
75 




(167r^)a A 
( 

TSn^ A 

yey-a 

{len^y^ A 

yaVuS^ 

(ISn-'y-i A 

ydyu 

(ISn-'y-' A 
„2 2 

yu V 
(le-iv'-'y-' A 

2 2 9 

via V-' 

(167r2)3 A 
y^/S (1671^)3 A 


( 1 4- ^ 
1 ,v^\ 

A^ y 

f 1 ,v^\ 
I 16^ A^ i 

^^ + 4^ 

f 1 , vA 
I 16^ ^ A^ y 

^^ + 4^ 

\Ta^ + A^j 




< 0.5 
2 X 10^ 
2 X 10^ 
7 

4 X 10^ 
2 X 10^ 
2 X 10^ 
1 


PPOv, mix, HElnv, ILC, LHC 

f3(30u, LHC 

PjSQv, mix 


D 

U 
U 

u 
u 
u 
u 
c 



The third column of Table [H labeled map, presents our estimate for the operator-induced Majorana neutrino mass 
expressions. These were derived based on the estimation procedure discussed earlier. Trivial order one factors, as well 
as the generation dependent coupling constants A have been omitted, as already advertised. Flavor specific charged 
lepton Yukawa couplings are explicitly denoted and to distinguish them from yi, y„ and yd, meant to represent 
a, /9-independent Yukawa couplings. A summation over all "internal flavors" is assumed for each entry. For order one 
coupling constants, this sum is strongly dominated by third generation Yukawa couplings. Upon setting these mass 
expressions equal to the observed scale of light neutrino masses (0.05 eV), we extract the required cutoff scale A for 
each operator. This quantity, defined to be Ai,, is listed in column four in units of one TeV. Numerical results were 
obtained assuming the current best fit values for all SM parameters. Associated errors are negligibly small as far as 
our aspirations are concerned. 

Fig. [5] displays the distribution of extracted cutoff scales, Aj^. The histogram bars are color coded to reflect the 
different operator mass dimensions. The distribution spans thirteen orders of magnitude, from the electroweak scale to 
10^^ TeV. It is interesting to note the general trend of operator dimension with scale: as expected, higher dimension 
operators are characterized by lower ultraviolet scales. For operators associated with the lowest ultraviolet cutoffs, the 
lepton number breaking physics occurs at the same energy scale as electroweak symmetry breaking. In this case, one 
needs to revisit some of the assumptions that go into obtaining the bounds and predictions discussed here. Regardless, 
it is fair to say that some of these effective operators should be severely constrained by other experimental probes, as 
will be discussed in the next section. 

The natural scale for most of the explored operators is well above 10 TeV, and thus outside the reach of future 
experimental efforts except, perhaps, those looking for neutrinoless double-beta decay. The remainder, however, 
should yield observable consequences in next-generation experiments. This small subset arguably contains the most 
interesting cases on purely economic grounds, as they naturally predict tiny neutrino masses as well as TeV scale new 
physics, which is already thought to exist for independent reasons. It is aesthetically pleasing to imagine that all, or 
at least most, of nature's current puzzles can arise from the same source, as opposed to postulating various solutions 
at different energy scales. It is important to note that one can "push" more of the operators into the observable TeV 
window by modifying the coupling of the new physics to different fcrmion generations. In particular, since many of 
the induced neutrino masses depend upon fermion Yukawa couplings, one can efficiently reduce scales by simply and 
uniformly decoupling the third generation. In most cases this can yield a Ki, reduction of several orders of magnitude; 
a factor that can be further enhanced by also decoupling the second generation. Under these conditions, the resulting 
distribution, analogous to Fig. [51 would show the majority of the operators piled up near and slightly above the 
electroweak scale. A detailed exploration of this possibility would be impractical and is not pursued further. We 
will, however, like to emphasize that this strategy of decoupling the new physics from the heavy fermions is very 
non-standard. In most cases, one is tempted to decouple light fermions from new physics both because these lead 
to the strongest constraints and because one tends to believe that the large Yukawa couplings of third generation 
fermions are entangled with the physics of electroweak symmetry breaking. 

Not all extracted cutoff scales are subject to a strong dependency on SM Yukawa couplings. In particular, the 
A^ values for the majority of dimension-eleven operators in the large histogram bar near 10^ TeV would not shift 
down at all under this hypothetical decoupling of the third generation from the new physics. These are the operators. 
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FIG. 2; A summary histogram of the scale of new physics extracted from the 129 LNV operators introduced in Table U 
We assume a radiatively generated neutrino mass of 0.05 eV and universal order one coupling constants. The contributions of 
operators of different mass dimensions are associated to different colors (shades of gray), as indicated in the caption. 



as shown in Diagram (d) of Fig. [l] whose induced neutrino mass matrix is independent of the Yukawa sector. In 
such cases, mct/3 are only functions of the various gauge couplings. As such, these constitute the most robust results 
of our analysis. These operators all predict an anarchic Majorana neutrino mass matrix of overall scale given by 
rriu = / [IQtt'^Yv^ / implying an energy scale ^ 10^ TeV. The only other "Yukawa invariant" cutoff scale 
estimate arises for the dimension- five operator Oi. Oi captures the physics of all versions of the seesaw mechanism [ol], 
and is at the heart of most of the model building currently done within the neutrino sector. Its ultraviolet completion 
can precede in only three distinct ways [13] . These possilailities are via the exchange of heavy gauge singlet fermions 
(type I seesaw), SU{2)l triplet scalars (type II seesaw) [ll|, SU{2)l triplet fermions (type III seesaw) 12 , or some 
combination thereof. Its popularity is well-founded for a number of reasons, including its underlying simplicity in 
structure as well as the purely empirical fact that it is the "lowest order means" of neutrino mass, and as such is 
easily generated by a "generic" LNV model. Additionally, the high scale associated with the seesaw mechanism can be 
easily incorporated within existing theoretical models and serves to help explain the observed baryon antisymmetry 
of the universe via leptogenesis [l3l| . 

For the purposes of direct observation, Oi's high cutoff scale, nearly 10^^ TeV, places it well outside of the "de- 
tectable region" (A < 10 TeV) and renders it uninteresting for the purposes of our analysis. Of course, there always 
remains the possibility that Oi is generate by very weakly coupled new physics (or very finely-tuned new physics 
[T^). in which case we expect to run into the new ultraviolet degrees of freedom at energies well below 10^^ TeV. In 
the case of Oi , it has been argued that new physics at almost any energy scale (from well below the sub-eV realm to 
well above the weak scale) will lead to light neutrino masses [lB| without contradicting current experimental results. 
Such possibilities - related to the fact that the new physics is very weakly coupled - are not being explored here, as 
we always assume that the new degrees of freedom are heavier than typical experimentally accessible energy scales. 

Armed with our derived new physics scales Ai/, we proceed to plug them back into the different irrelevant LNV 
operators and search for possible means of future observation as well as already existing constraints. Generally, those 
operators that yield the largest experimental signals have the lowest cutoff scales. We conclude that, if associated to 
neutrino masses, the effective cutoff scale A^ of the following effective operators is constrained to be less than 1 TeV: 



C34, C36, 037, 03S, C53, 054a^b,c,d ' ^56, C59, Ceo, C'to- 



(II.7) 
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These may lead to observable effects at future high energy accelerator facilities. Additionally, such low scales may 
also indirectly lead to observable effects in "low energy" (but high sensitivity) experiments. There are more operators 
associated with slightly higher scales between (1 — 10) TeV that may manifest themselves experimentally via virtual 
effects. These are 

OiQ,On, Oi8, ©19, O35, ©43^^^ , ©50, ©52, ©SSa.b.c: ©58, ©65, ©68„.6 , ©73a, ©75- (H-^) 

These operators yield finite predictions for more than one observable, such that experimental efforts in seemingly 
unrelated fields can help constrain the class of possible LNV models or even help identify the true LNV model. 

III. GENERAL OPERATOR CONSTRAINTS AND PREDICTIONS 

There are, currently, bounds on LNV processes from a number of independent experimental sources (igI. [l7j. Many 
of these are presently too mild to constrain the operators listed in Table U once their ultraviolet cutoffs A are set to 
the required value indicated by the presence of non-zero neutrino masses, A^. The situation, however, is expected 
to improve in the next several years with increased rare decay sensitivities and higher collider energies. Here we 
survey the experimental signatures of these operators in terms of the minimal scenarios described above. Specifically, 
we address the potential of neutrinoless double-beta decay (Sec. IIII Ap . rare meson decays (Sec. IIIIB[) . and collider 
experiments (Sec. IIII Cp to constrain the effective operators in question, assuming that, indeed, they are responsible 
for the observed non-zero neutrino masses. As before, we will use the approximations discussed in Sec. |TT1 and warn 
readers that all the results presented are to be understood as order of magnitude estimates. The results, however, are 
useful as far as recognizing the most promising LNV probes and identifying different scenarios that may be probed 
by combinations of different LNV searches. 

Most of this section will be devoted to probes of LNV via simple variants of the following process, which can be 
written schematically as 

ialp ^ d^d^UpU^. (III.l) 

Greek subscripts run over all different fermion flavors. Given the assumed democratic models, coupled with our present 
lack of experimental information, one would expect that all flavor combinations are equivalent to zeroth order. Any 
indication to the contrary would signify important deviations from simple expectations, and thus begin to reveal 
the flavor structure of the new physics. The above selected "golden modes" often yield the largest LNV rates, but 
this is not always the case. For example, some operators do not allow tree-level charged dilepton events, but rather 
prefer to include neutrino initial or final states. LNV processes with initial and final state neutrinos are extremely 
difficult to identify. The only hope of such discovery channels is, perhaps, via neutrino scattering experiments on either 
electron or nucleon targets, using well understood neutrino beams. We point out that any neutrino/anti-neutrino cross 
contamination induces ambiguity onto the total lepton number of the incident beam and would serve as a crippling 
source of background for LNV searches. This reasoning rules out conventional superbeam fTs*! facilities as well as 
proposed neutrino factories |.19i] . which contain both neutrino and anti-neutrino components, but does suggest modest 
possibilities for future beta-beams (20| . Given projected beta-beam luminosities and energies along with the derived 
cutoff scales A^,, it seems unlikely that LNV can be observed in such experiments. Another possible discovery mode 
involves only two external state quarks and an associated gauge boson as in the sample process £a(f3 du + W~ . It 
turns out that the rates for such processes are generally suppressed for the majority of operators involving six fermion 
fields, as we are trading a phase space suppression for a stronger loop suppression. For those operators with only four 
fermion fields, the situation is not as straightforward and, in some cases, the three particle final state is preferred. 
Typically, the neutrino mass induced cutoff scales of those operators are high {A^ » 100 TeV), so it would be quite 
difficult to observe such effects. Of course, any ly-boson final state will either promptly decay leptonically, yielding 
missing energy and unknown total lepton number, or hadronically, reducing the reaction back to that of the golden 
mode. 

Another possibility is to replace two or more of the external quark states in Eq. pil.ip with leptons in such a 
way as to preserve charge, baryon number, and AL = 2 constraints. While many operators favor this structure, a 
little thought reveals that at least one external neutrino state is always present, which leaves only a missing energy 
signature, and little means of lepton number identification in a detector. Such events would not be clean, but of 
course, three final-state charged leptons and missing energy are enough to extract the existence of at least AL = 1 
LNV, provided that the number of invisible states is known to be no greater than one. This last requirement is difficult 
to achieve in the presence of the large backgrounds and the limited statistics expected at future collider facilities, 
but should still be possible given a concrete model probed near resonance (see for example '21']). Therefore, while 
important and potentially observable, this mode is not generally the best place to look for LNV and is neglected in 
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the remainder of our analysis. From this perspective, the only other relevant channel of LNV discovery is related 
to W and Z rare decays into final states with non-zero total lepton number. This possibility is briefly addressed in 
SeciniBl 

A. Neutrinoless Double-Beta Decay 

Here we probe the expectations for neutrinoless double-beta decay {f3(3Qi^) for each operator listed in Table U (3/30iy 
is the LNV (AL ~ 2) process where, within a nucleus, two down quarks convert into two up quarks with the emission 
of two electrons but no neutrinos, or in the language of nuclear physics {A,Z) {A,Z -I- 2) -I- e~e~. See [l^l and 
references therein for a comprehensive review. While precise computations of nuclear matrix elements are essential 
for making detailed predictions [l^l , the minimal parton-level description given above is adequate for the purposes of 
this study. There is a continuing legacy of cutting edge experiments designed to search for (jpOv with no success to 
date.* Currently the ^^Ge half-life for this process is bounded to be greater than 1.9 x 10^^ yr and 1.57 x 10^^ yr 
at 90% confidence by the Heidelberg-Moscow f25| and IGEX [2()| experiments, respectively. Future experiments are 
poised to improve these limits (for several different nuclei) by a couple of orders of magnitude within the next five to 
ten years (27| . 

If one assumes that (3(30^ proceeds via the exchange of light Majorana neutrinos, its amplitude is proportional to 
the ee element of the Majorana neutrino mass matrix, 

3 

mee=^m,f/2, (IIL2) 

4=1 

where rrii are the neutrino masses and Uei are elements of the leptonic mixing matrix. With this, one can extract 
the upper bound mee < 0.35 eV (90% confidence level bound, 17]) from current experiments while next-generation 
experiments are aiming at niee ^ 0.05 eV^ [131 . In general, LNV new physics will lead to additional contributions to 
fipOv, most of which are not proportional to ruee- However, the amplitude for (3(301/ can still be expressed in terms of 
an effective mee, m-gf , which is an operator-specific quantity that will be used to analyze new models of LNV. 

Here, we define six different "classes" of diagrams one can construct out of LNV irrelevant operators that contribute 
to (3(3Qv at the parton level. These are illustrated in Fig. [31 and classified by the dimension of the generated LNV 
interaction, depicted by large gray dots. In order to unambiguously separate the different classes, note that the grey 
circles are defined in such a way that all fermion and Higgs legs that come out of it are part of the "parent" operator 
O (and not attached on via reducible SM vertices), while all other interactions are SM vertices. The dots should be 
viewed as hiding the underlying LNV interactions. In general, they contain a mixture of coupling constants and loop 
factors that must be evaluated explicitly for each diagram. It is important to emphasize that the contribution of a 
generic operator O to (3(30v will consist of contributions from all different classes, while usually dominated by one of 
them. We show the lepton number conserving electroweak vertices (point-like) as effective four- fermion interactions, 
justified by the low energy scale of nuclear beta decays. The dotted lines indicate the exchange of VF-bosons, labeled 
by W and H (charged Higgs goldstone boson). Helicity arrows are explicitly included where uniquely determined, 
implying that the arrowless legs can have any helicity. 

describes the standard scenario of (3fiQv mediated by light Majorana neutrinos. It is simply two electroweak 
vertices held together by a Majorana mass term on which two neutrinos are annihilated. The amplitude for this 
diagram is proportional to ra^e, as defined in Eq. pil.2[) . The dependence on such a neutrino mass is intuitively 
clear considering the need for a helicity flip on the internal neutrino line. The remaining diagrams are qualitatively 
different from this standard case. Most importantly, none of them require "helicity flips" and are therefore not 
directly proportional to neutrino masses. They are, however, proportional to inverse powers of the new mass scale 
A^, and hence also suppressed. These effects are not entirely independent, since the value of A^ was extracted from 
the requirement that neutrino masses are small, but correlations are relaxed enough to allow nontrivial consequences. 
It is this partial decoupling from neutrino masses that allows larger than naively expected contributions to (3(3Qv 
from some of the LNV irrelevant operators. Before proceeding, we make the trivial observation that the amplitudes 
following from D^,, D4 and 1)5 are additionally proportional to two powers of CKM matrix elements, namely iKidP, 



* There is currently a positive report of /Bf^Ou at the 4.2(7 level by a subset of the Heidelberg-Moscow collaboration |24j . They report a 
measured half-life of 1.74l^;i^ X 10^^ years which maps to ~ (0.2 — 0.6) eV. We choose to neglect this controversial result, which 
is still awaiting independent conformation. 

t The parameter change from half-life to mee depends heavily on nuclear matrix element calculations. Current calculations induce an 
uncertainty of less then a factor of four on mee for most parent isotopes |23|| . 
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FIG. 3: The parton level Feynman diagrams contributing to neutrinoless double-beta decay, labeled by the dimension of the 
underlying lepton number violating interaction, indicated by gray dots. Each diagram is generated, at some order in perturbation 
theory, by all analyzed interactions, but estimates of their magnitudes depend heavily on the details of the operators, including 
their associated scale A^, fermion content and helicity structure. 



whereas Dq and Dj are only proportional to one power of V^d-^ The tree-level diagram Dg has no CKM "suppression." 
While this is a purely academic fact in the case of Pf30iy {\Vud\ ~ 1!)) it leads to important consequences for analogous 
rare decays that depend on the much smaller off-diagonal CKM matrix elements. We will return to these in the next 
subsection. 

For a given operator, the relative size of each diagram's contribution to the total decay rate depends on many 
factors including the operator's dimension, scale, fermion content and helicity structure. The dominant contributions 
must be calculated on a case by case basis. Generally, the high scale operators (A^ > 10 TeV) are dominated by 
the two dimension-four diagrams and since many factors of A,^ will be canceled by divergent loops inside the 
gray dots thereby minimizing the suppression. All else being equal, D^, is the strongest of the pair since it is 

enhanced by ~ from the two propagating neutrino lines as opposed to only ~ Q^^ for the one neutrino case 
shown in diagram D4. For those operators with no tree- level field content, can still be very important, but its 
dominance is nevertheless rare. As discussed in Sec. [Ill these are precisely the operators that have the greatest loop 
suppressions and consequently lower energy scales suggesting the need for diagrams beyond D4. The effects of low 
cutoff scale operators (A^ < 1 TeV) are not severely suppressed by l/A^, (by definition), so the dominant diagrams 
will typically be of the highest dimension allowed by the tree-level structure of the operator. For such low scales 
and for operators of the following schematic form dduuee (dimension 9) or dduueeHoHo (dimension 11), Dg always 
dominates the /3/30i^ rate yielding amplitudes proportional to 1/A^ and w^/A^, respectively. For intermediate scales, 
and when the operator's field content does not directly support P(30iy due to lack of quark fields, the situation is 
not as straight forward and one must perform the relevant computations to determine the dominant diagrams. Still, 
it should be noted that diagrams containing internally propagating neutrinos are enhanced by inverse powers of Q 
and maintain a slight advantage over their neutrinoless counterparts. One can thus generally expect diagram Dq to 
dominate the decay rates for low scale operators when Dg is suppressed. The opposite is true for interactions 
taking place at higher energies in, for example, next-generation colliders, as discussed in Sec. IIII CI 

Since each diagram in Fig. [3] can have different external helicity structures, the different contributions to the total 
rate will be added incoherently, thus eliminating the effects of interference. There are some case specific coherent 
contributions that we neglected in our treatment since most rates are dominated by a single diagram. Another potential 
difference among the different contributions is related to nuclear matrix element calculations: can the calculations 
done assuming (3(30^ via the standard light Majorana neutrino exchange scenario of diagram Dy be applied to the 
more general cases encountered here? We have nothing to add to this discussion except to naively note that there is 
no obvious reason why such rates should be severely suppressed or enhanced relative to the standard scenario. We 
therefore assume that all nuclear matrix elements are identical and can be factored out of the incoherent sum. We 



t This is true provided that we assume no flavor structure for the underlying operator, or, equivalently, that all dimensionless coupling 
constants are order one. If one is motivated by experiment to postulate a minimally flavor violating scenario, to perhaps ease constraints 
from flavor changing neutral currents, the statement must be modified accordingly. 
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assume that this approximation is not more uncertain than the other sources of uncertainty inherent to our study 
(Ukely a very safe assumption). 

As drawn, each diagram Di contributes to the amplitude that characterizes j3[3Qv. For example, the amplitude 
associated with D^j is proportional to 



where Gp is the Fermi constant. The remaining diagrams will contribute with Aoi oc C,{v, up to a dimension- 

less coefficient containing various numerical/loop factors, as well as general scale depencies parameterized by some 
power of the ratio w/A. The function C,(v, Q) has mass dimension i — 9 so that all Aoi have the same mass dimension. 
Note that all aspects of Aoi are calculable given a LNV operator and diagram. We can analyze each operator in 
terms of an effective rrf^, defined in terms of the underlying dimension nine amplitude Aoi by 



where i runs over the set {i/, 4, 5, 6, 7, 9} that labels the diagrams shown in Fig. [3l and Q ~ 50 MeV is the typical 
momentum transfer in (3(301/. m^f can be directly compared with experiment and used to make prediction for future 
observations. A few comments are in order regarding this quantity. First, it is a useful derived object that has no 
direct connection to a real neutrino mass and is valid to arbitrarily large values. Note that in the case of Majorana 
neutrino exchange, rrigf = rriee only if iriee ^ Q- When neutrino masses are greater than Q, m^f cx 1/m. Our 
definition of m^f also conforms to the use of large effective masses in [l^. The second comment is that, unlike the 
case of niee, which is valid for any process involving the exchange of electron-like Majorana neutrinos, is case 
specific. It must be calculated separately for each process, as each one, in general, is composed of different diagrams. 
In particular, the calculations of the effective mass for (3(301/ expressed here are not directly applicable to other LNV 
processes and should not be interpreted as such. 

The rrigg distribution extracted from all operators is shown in Fig. [?] assuming the scales A^, derived in Sec. |TT] 
and color-coded for convenience within the histogram. Specifically, we indicate in green the operators that are 
characterized by sub-TeV scales and thus accessible to next-generation experiments via direct production. The blue 
and red operators are characterized by scales between (1 — 5) TeV and (5 — 25) TeV respectively, where virtual effects 
should be most important for collider searches. The majority of operators, shown in cyan, are suppressed by scales 
greater than 25 TeV and are hence quite difficult to observe in other search modes. We also explicitly label each 
operator within the histogram bars for easy identification and comparison. One should notice the expected general 
trend that increasing leads to a decrease in m^g and vice-versa. The vertical axis is truncated at 15 operators, as 
the bar near 0.05 eV, dominated by the light Majorana neutrino exchange described above, would extend to nearly 
100 operators. With broken vertical lines, we indicate the current 90% upper bound [l3|, "T-ef = 0-35 eV, and the 
potential reach of future experiments. 

This distribution, which spans over six orders of magnitude (from 10^"* eV till 10^ eV), reveals many important 
features of the effective operator set. Beginning at the largest values, we find that the twelve operators appearing 
near 300 eV all have the expected common feature of low energy scales, including Oig with A^, only just above the 
1 TeV mark. Additionally, the contribution of the majority of these operators to (3(301/ is dominated by the tree-level 
Dg diagram. The exceptions are ©54^ ^ and ©70, all of which are characterized by sub 0.5 TeV scales and dominated 
by diagram Dq. Consequently, these are subject to a loop and Yukawa/gauge^ suppression relative to their Dg 
dominated cousins, but the difference is not visible given the resolution of the figure. It is interesting to note that 
these three operators have the correct quark and lepton content for large (3(30v, but their SU{2)l gauge structures 
forbid large tree-level contributions. Similarly, operators Oie, O^^^ (, ^, CeSaJ and ©75 are also dominated by diagram 
Dq accompanied by slightly higher cutoff scales. This drives down m^^ significantly considering the leading one-loop 
scale suppressions of A^^ and A^'' for the dimension-eleven and dimension-nine operators respectively. We point out 
that operators ©54^ ^ ^ ^ and ©55^ ^ ^ yield almost identical expressions for their respective (3(30v amplitudes (as well as 
their radiatively generated neutrino mass expressions) with up and down quark Yukawa couplings exchanged. While 
this action enhances most of the ©55 (3(30v couplings relative to those of O54, it also raises the 055 A^, scale by nearly 
a factor of four and thus drives mgf down by orders of magnitude. 



\Vud 




(IIL3) 




(IIL4) 



5i As it turns out these are all suppressed by a single bottom quark Yukawa coupling as well as two powers of the SU(2)]^ gauge coupling 
g, but this fact cannot be deduced from Fig. |4] alone. 
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FIG. 4: TTigf distribution derived for the neutrinoless double-beta decay process as described in the text. The calculations were 
made assuming the scales derived in Sec. [Ill as well as universally order one coupling constants. The histogram bars are 
labeled explicitly with operator names and color-coded by their cutoff scales. Also shown in light gray is the region probed by 
next-generation experiments. The vertical axis is truncated at 15 operators to best display the relevant features of the plot. 



The remaining operators all predict < 1 eV, close to current experimental bounds. The histogram bar near 
0.1 eV is composed of operators of very different scales. ©34 and ©56 are both characterized by low cutoff energy 
scales around 0.5 TeV, but, due to their fermion and helicity structure, their contributions to P(30iy are dominated by 
two- loop versions of diagram Dq. The neutrino-mass-required cutoff for ©73^ is around 7 TeV and its contribution 
to f^pOv is also dominated by diagram Dq. In this case, however, the two-loop version turns out to be larger than 
the allowed one-loop amplitude due to strong scale suppressions (the added loop reduces the cutoff dependency from 
A^''' to A^'^). This behavior is characteristic of operators with a larger value of A^. The A^ = 40 TeV operator ©20, 
defines the lower edge of this histogram bar. It is dominated by the one-loop diagram Dq enhanced by a top quark 
Yukawa coupling and, being a dimension nine operator, is only suppressed by A^^ from the start. The next bar down 
contains operators dominated by D^. Most of these are suppressed by a very high energy scale, but a small subset 
is characterized by scales A^, < 25 TeV. In particular operators On, Ois and ©57 are all cutoff at 2 TeV but, due 
to their fermion content they cannot participate in any of the non-standard interactions of Fig. [3] at a low enough 
order in perturbation theory. Similarly, the intermediately scaled operators 043^^^, O50, O52, 062, Oesi and Ogg^ 
have either the wrong fermion content or gauge structure to enhance any of the (3(30^ diagrams (other than D^,) to 
an observable level. These operators are important because their minimal forms are experimentally unconstrained 
yet still potentially observable to both next-generation j3f3Qv and collider experiments. The remaining histogram bars 
with m^f < 10-2 eV are not accessible to /3/30i^ experiments in the foreseeable future. Each of these diagrams are 
dominated by D^, either due to high suppression scales as in the case of O7 and Ogj or, as in O35 and O58, the 
operator's fermion content simply disfavors other contributions to the j3[3Qv amplitude. It is the general form of the 
neutrino mass matrix derived in Table IH where we see that niee oc ye, that drives these operators away from their 
peers near m^f = 0.05 eV. It is unfortunate that the two "low" dimensionality operators O7 and Og are cutoff by 
energy scales in excess of 100 TeV and are hence invisible to any direct probe. If either of these operators have 
anything to do with nature, it is unlikely that LNV will be observed in the foreseeable future in any experiment. 
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On the other hand, any observation of LNV will rule out these types of scenarios. Additionally, as will become clear 
shortly in Sec. IIV[ current neutrino oscillation data already marginally disfavor such operators and have ample room 
to tighten constraints in the near future. 

It is interesting to point out that the lower boundary of the currently excluded region falls within the 771°^ dis- 
tribution, suggesting exciting prospects for the future. That being said, one should not read too much into current 
and future null results as, for most operators, relatively small cancelations and order one factors, not accounted for 
here, can push the relevant rates below the observable level depending on the underlying ultraviolet theory. On the 
other hand, one is allowed to interpret that operators that lead to nfj^ > 10 eV are severely constrained (if not ruled 
out) as proper explanations for neutrino masses if one assumes the new physics to be flavor "indifferent" - order 
one factors cannot be evoked to save the scenario. Once this assumption is dropped, however, it is quite easy to 
"fix" these scenarios, since the large /9/30i/ rate is a direct consequence of the universal order one couplings and the 
relatively low cutoff energy scale A^. For example, one can suppress the coupling of new physics to first generation 
fermions (compared to second and third generation fermions), thereby suppressing the worrisome diagrams of Fig.[3l 
This will have little effect on the relation between A and the neutrino masses, discussed in Sec. HIl since these are 
either generation independent or highly reliant on third generation Yukawa couplings. Of course, by combining fSpOv 
searches with other probes we can obtain a much better idea of the origins of LNV as well as the relevant model(s), 
if any, chosen by nature. 

B. Other Rare LNV Decay Processes 

Most of the qualitative discussions of Sec. IIII A| devoted to (3f30i', can be directly applied to other rare decay 
processes with the same underlying kernel interaction described by Eq. pil.ip . For such processes one need only 
analyze simple variants of the diagrams listed in Fig. [31 using crossing amplitude symmetries to account for the 
needed initial and final state fermions. Other factors must be added to the various electroweak vertices to account for 
quark flavor mixing. The requisite CKM matrix elements can highly suppress many diagrams for processes involving 
cross-generational quark couplings. In fact, only tree-level Dg diagrams are safe from such suppressions. Next, and 
most importantly, one must include the appropriate characteristic momentum transfer Q of the new system. Specific 
rates are highly dependent on this quantity as effective operator cross-sections typically grow with some power of 
Q. The particular exponent of the power law depends on the diagram, but naive dimensional analysis dictates that 
r cx Q^^ for diagram Dg, rendering it highly dependent on a reaction's energy transfer. The fact that each diagram 
varies with Q in a different way implies that predicting the dominant contributions to a given process is non-trivial 
and must be addressed quantitatively. Finally, in the cases of hadronic decays, one must also account for initial/final 
state matrix elements. We assume that all factors can be simply estimated on dimensional grounds. 

Unlike the f3/30i' case, some meson decay modes proceeding via new LNV tensor interactions are expected to be 
suppressed. Such processes are one instance in our analysis where an operator's Lorentz structure can qualitatively 
affect expected LNV decay rates. One can understand this by considering a meson decay mediated by a new tensor 
particle. The parton level interaction has the form (ua^^c?)T'"^ where the initial state quarks are explicitly shown 
and all other fields are contained in the tensor T'^" . Following the standard procedure we factor out the hadronic 
structure in the form of a free decay constant and write the amplitude as generally allowed by Lorentz invariance in 
terms of the external state's four-momentum. Due to the antisymmetry of u^i/, this amplitude vanishes to first order. 
Non-zero contributions to this decay mode must necessarily involve individual parton momenta and are therefore 
suppressed relative to the usual vector-like decay calculations. From this, it is clear that models of LNV containing 
tensor couplings will often evade the predictions and bounds of this section. Tensor operators will mediate LNV meson 
decays into more complicated final states (one may include, say, initial/final state radiation). Associated rates are, 
however, subject to additional gauge coupling and phase space suppression that tend to further reduce the already 
tiny LNV rates beyond any hope of detection. 

Rare LNV meson decays have been experimentally pursued for many years (171] . Here, we focus on the AL = 2 
processes M' ^ M + t^P^-, where M' and M are the initial and final states mesons respectively and the Is represent 
like-sign lepton pairs of arbitrary flavor. Electric charge conservation dictates that M' and M have equal and opposite 
charge. Here we take each meson to consist of a color singlet up-type/antidown-type bound state* and factor out all 
long distance hadronic effects. In this way we can view the meson decay process as du — > ^af/j -I- du for all up-type 
and down-type quark flavor combinations. The effective LNV diagrams contributing to this process are shown in Fig. 



* For simplicity we assume that both the process M' — > M -|- iaf-ji and its conjugate have similar amplitudes and therefore treat them 
symmetrically. Large CP-violating effects can invalidate this assumption. 
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FIG. 5: The parton level Feynman diagrams contributing to rare LNV meson decay labeled by the dimension of the underlying 
lepton number violating interaction, indicated by gray dots. Each diagram is generated, at some order in perturbation theory, 
by all analyzed interaction, but estimates of their magnitudes depend heavily on the details of the operators, including their 
associated scale A^, fermion content, and helicity structure. 



[5] with the same naming scheme as their analogs in Fig.[3l Here, V and V' denote potentially distinct elements of the 
quark mixing matrix. We additionally point out the potential dependency on all entries of the Majorana neutrino 
mass matrix elements rriafs in diagram D^, as opposed to the PfiOv case where Dn depends only on ruee- These 
processes probe combinations of the neutrino masses that are naively unconstrained by /3/30;/ [2^. In general, the 
varied flavor structures encountered in meson decays allow for experimental probes into new physics couplings across 
the fermion generations. We pointed out earlier that some of the LNV operators lead to unacceptably large rates for 
PfjQv unless first generation quarks participate in the new interactions with severely suppressed couplings (compared 
with second and third generation quarks). If such a scenario is realized in nature, rare D or B decays may be much 
more frequent than naive expectations. For this reason, improving rare decay sensitivities to all channels is essential 
to completely constrain models of new LNV physics beyond the minimal framework analyzed here. 

Reference [l6j summarizes LNV upper bounds on all of these processes in terms of the effective Majorana neutrino 
mass matrix element rn°^ that one would extract from observation assuming that all decay rates are dominated by 
the light neutrino exchange shown in D^. Hence, we can compare operator expectations with current experimental 
limits in exactly the same way as was done in Sec. IIII Al For a given LNV meson decay, m°*g is defined from the 
contribution of the different classes of diagrams to the rare meson decay in question, exactly as rrf^ was defined in 
the previous subsection (see Eqs. piI.3|III.4|) ). Direct estimates for different process reveal to°^ distributions similar 
to that for m^f depicted Fig. lU up to "rescalings" that reflect the different kinematics and the presence of small 
CKM mixing matrix elements. Results are summarized in Fig.[n]for a representative sample of charged meson decays. 
Each histogram is labeled by its associated decay mode and is color-coded to indicate the neutrino-mass constrained 
cutoff scale of the different LNV effective operators. For simplicity, we refrain from listing operator names on 
the individual histogram bars (as opposed to what was done in Fig. |4|). The "operator ordering" is very similar to 
that of Fig. 21 especially in the low A^ scale, high effective mass regime where decay rate predictions are particularly 
important. Note that the horizontal axes are relatively fixed for easy comparison and that the vertical direction is 
truncated and does not reflect the true "height" of the lowest mass bar (order one hundred operators). 

Specifically, we present effective Majorana neutrino mass distributions for the processes, reading down the panels 
from left to right, D ^ t: + D^K + €±£±, iP, ^ tt + ^±4, D^^K^ £±4, B ^ tt + B^K + ^±4, 

K 7r + ^^^^, as well as the rare r decay, M M' -\- J Here the final state leptons can be of any flavor allowed 

by energy conservation. Since, as previously discussed and explicitly verified numerically, the specific details of the 
distributions are mainly dictated by kinematics and CKM matrix elements, these results are robust under changes in 
the final state lepton flavors. The r decay distribution shown in the lower right panel is representative of all possible 
decay products including flrst and second generation charged leptons and light meson states. One should notice 
the expected general operator trend within each histogram as the characteristic cutoff scale is decreased, as well as 



t The actual calculations displayed in Fig. |6] assumed the charge lepton flavors la^p = Aie, while the r decay histogram (lower right-hand 
panel) was produced assuming the final state mesons MM' = KK. 
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FIG. 6: m'^p distribution for several rare LNV meson and r decays. Calculations assumed the charge lepton flavors £a£fi ~ f^e, 
while the r decay histogram (lower right-hand panel) was obtained assuming the final state mesons MM' = KK. The histogram 
bars are color-coded by suppression scale. Current bounds on these processes are typically above 1 TeV and are not visible at 
these small scales. 



the expected peaks near 0.05 eV dominated by light Majorana neutrino exchange. Additionally, each distribution is 
much "broader" than the one in Fig. 01 This observation exemplifies the fact that effective mass calculations depend 
critically on the underlying process. Indeed, maximum m^f values can reach nearly 10^° eV for the — > K~ +e^e^ 
decay but only lO'^ eV for PPQv. Current upper bounds for m'^^ from these processes, mostly well above one TeV, 
are well beyond the largest operator predictions here, ranging from mf^ < 0.09 TeV for the case of if+ n e+/x+ 
to < 1800 TeV for the case of K~ fi^ [l^. It is curious that the best meson decay bounds come from 

rare LNV kaon process but, as can be seen in the lower left panel of Fig. [6l these yield by far the lowest predictions. 
Future experiments have the potential for observing LNV for a select few operators only provided vast improvements 
in meson production luminosities. Current and upgraded B-factories [l^ are expected to provide the most significant 
improvements, considering the large derived i?-meson effective masses shown in Fig. [HI Still, the best cases from the 
figure yield only the tiny branching fraction 1.8 x 10~^^ for the case of the rare decay B'^ — s- 7r~e"'"/i"'", nearly eleven 
orders of magnitude below the current experimental limit of 1.3 x 10~^ 

Another possible search mode involves the decay of the Z-boson into LNV final states. The dominant contributions 
to this process are generally unrelated to the reactions summarized in Eq. (jIII.l[) and shown schematically in Figs. [3] 
and m While there is a slight connection between them as one can always attach a Z-boson to various fermion lines 
in each diagram, there are potentially large lower order contributions arising within the operators themselves. The 
latter, when present, can easily overtake the associated "golden mode" counterparts. In this context, such processes 
can be thought of as the decay of the longitudinally polarized Z-boson. Strict bounds exist on such decays from 
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the LEP-I [sOl and SLD 'si'l experiments. Each element of the operator set predicts decays into final state fermions 
with total lepton number L — 2. The dilepton pair can be of any flavor and is generally accompanied by two or four 
additional fermion states, depending on the dimension of the operator. We restrict our discussion to the dimension- 
eleven operators comprising the majority of the sample, as these are typically suppressed by lower cutoff A^, scales 
and, equally important, explicitly contain Higgs doublets in their field content. In this case, tree-level decays result 
in a six-fermion final state which suffers from a large phase space suppression and cumbersome multiplicities that are 
likely to render even the most sophisticated search ineffective. The only possibility of this type that yields a charged 
dilepton signal is Z i^t^qqqq (quarks of all allowed flavors implied), but many other possibilities exist involving 
invisible final state neutrinos. A little thought also reveals that closing fermion loops in an attempt to obtain simpler 
final states and thus render the analysis more tractable will necessarily result in final state neutrinos. Therefore, 
the majority of the Z-boson LNV decay channels involve invisible final states with practically undetermined total 
lepton number. The prospect of direct discovery by these means seems dismal, but indirect constraints on LNV are 
still possible from bounds on the Z-boson invisible decay width. There is currently a statistically insignificant, but 
nonetheless captivatin g, 2 (j deviation between the observed invisible decay width and SM expectations assuming three 
light neutrino species [I7l|. The experimentally extracted branching ratio was found to be slightly smaller than its 
predicted value so that a new LNV contribution of the form Z —f Vai^i3 would push the invisible branching ratio 
in the "wrong" direction. From these bounds the decay width of any new contribution to the Z-boson decay is 
constrained to be less than 2.0 MeV at the 95% confidence level [13) 113 ■ ^ quick estimate reveals that this constrains 
the dominant LNV amplitudes Az < \/47r(2.0 MeV/M^) ~ 0.53. For the dimension-eleven operators of interest, the 
largest possible amplitude is of order ?/^/(I67r^)^(w/A)^ where y is an arbitrary fermion Yukawa coupling and four 
powers of the cutoff scale A are removed by divergences in the closed diagram loops. The constraint above translates 
into y'^{v/A)^ < 4.1 x 10^, which is easily evaded by even the best case scenario of y = j/t « 1 and A « u. Experimental 
bounds on r.m„ must be improved by a factor of a million before they start significantly constraining LNV (under the 
assumptions made here). This result holds for virtually all possible flavor structures. We conclude that rare Z-boson 
decays are not practical discovery modes for the LNV effects considered here, but look to future rare Z-boson decay 
studies for more information. 

In a similar way, one can also dismiss the case of rare W^-boson decays as promising probes of LNV. As in the 
Z-boson case, the W^-boson can decay into a variety of L = 2 final states proceeding either through couplings to 
left-handed fermion lines or explicit operator content. Here, however, there is no six-fermion, same-sign dilepton final 
state with no neutrinos due to conservation of charge and weak isospin, so the lowest order observable mode is already 
loop suppressed to W" ^ ^^^^ + qq. Current W^-boson decay bounds are far too weak to constrain such suppressed 
LNV T?| and are not likely to improve to the level implied by the operators under consideration, which predict the 
tiny decay rate Tlnv < mn/(47r)/(167r^)'^(u/A)^'^ « lO"*^ MeV in the best case scenario of electroweak scale A^. 
We also point out that, contrary to the Z-boson decay limits, there are no robust, indirect bounds that can be used 
to constrain LNV in the case of the VF-boson. Note that, despite dismal prospects for gauge boson decay driven 
LNV discovery within the minimal framework of "natural" effective operators, one can still construct theoretically 
well-motivated models that will yield observable signals. Particularly, in a weak-scale seesaw mechanism (Oi), the 
new degrees of freedom, comprised mostly of Majorana gauge singlet fermions (right-handed neutrinos), can mediate 
visible, AL = 2, W-boson mediated processes with little or no scale/loop suppression. This class of model is analyzed 
in (ssj and is exempt from the discussion outlined here. 

C. Collider LNV Signatures 

If neutrino masses are a consequence of ultraviolet physics related to cutoff scales around the TeV scale, we expect 
future high energy collider searches to directly access the new LNV physics. For example, the direct, resonant, 
production of new states could lead to rather spectacular signals of these models. It would also indicate the breakdown 
of the effective field theory approach undertaken here. To pursue such possibilities, one must assume a specific 
ultraviolet sector and study its signatures and implications on a case by case basis. In the looming shadow of the 
LHC, {34| and the more distant ILC [35], such an analysis is highly warranted but will not be pursued here. Instead, 
we assume that the masses of new ultraviolet degrees of freedom remain out of the reach of next-generation accelerator 
experiments. Such a situation can be easily accommodated within the context of the preceding results, considering 
the order of magnitude nature of the A^, estimates. 

We will concentrate on the process e~e~ — > qqqq (which will usually manifest themselves as jets) with no missing 
energy in an ILC-like environment [sH] with a center-of-mass energy of 1 TeV and an integrated luminosity of 100 fb^^. 
We also make the oversimplifying assumption that the detector system has equal acceptance to all quark flavors, and 
the ability to efficiently distinguish quarks, gluons and rs. By summing over all possible quark final states it is simple 
to estimate the total LNV cross section for each effective operator, assuming it is responsible for neutrino masses. Such 
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searches can be complemented by looking at e~e~ W~W~ , which have been discussed in detail in the literature 
[3^ . As discussed in Sec. IIIIBi the different LNV operators couple to one or more gauge bosons via an appropriately 
closed fermion loop or direct coupling to the Higgs doublet field. 

Charge and baryon number conservation dictate that the two quarks in e~e~ —^ qqqq are down-type quarks, while 
the two antiquarks are up-type antiquarks. At the parton level, the scattering process is similar to (3(30^, which 
motivates exploiting simple variations of the diagrams in Fig. [3] in order to calculate the relevant amplitudes, as was 
done in Sec. IIII Bl Here, the extensions are obvious: use crossing symmetry to rotate all lepton lines into the initial 
state and all quark lines to the final state taking special care to insert appropriate CKM matrix elements where 
needed. Due to the large characteristic momentum transfer Q of the e~e~ scattering, one must also "expand" the 
electroweak vertices and account for gauge boson propagation. With this in mind, the amplitude calculations can 
be carried over directly from the previous sections. Specific results are, however, quite distinct due to the higher 
center-of-mass energies involved. In the language of the underlying diagrams mediating this reaction, for diagrams 
characterized by TeV cutoff scales, diagram Dg, if allowed at tree-level, will dominate the rates. As in the previous 
cases, for intermediate to high cutoff scales, general diagram dominance must be addressed on a case by case basis. 
It is important to appreciate that, since these are non-renormalizable effective interactions, cross-sections grow with 
center-of-mass energy. For this reason, we expect many of the low cutoff scale operators to yield observably large 
signals at the ILC. 

Fig. [7] shows the e^e^ ^ qqqq cross-section distribution, in femtobarns, at the ILC, calculated for all 129 of the 
analyzed LNV operators. Once again, the extracted value of the cutoff energy scale A^ assuming constraints from 
neutrino masses are color-coded to indicate operators associated with a low (A,y < 10 TeV) or high (Aj, > 10 TeV) 
ultraviolet cutoff. Each bar is also labeled with the respective constituent operators, for convenience. Note that 
the vertical axis is truncated at fifteen operators (the left-most bin is over 60 operators high) to help clearly display 
relevant features of the plot. We also highlight the potential reach (defined as cross-section greater than the inverse 
of the integrated luminosity) of the ILC with a broken vertical line, assuming 100 fb~^ of integrated luminosity. 
This particular ILC luminosity value should be considered as a loose lower bound, introduced to give a feeling for 
the observable scales involved. It has recently been argued, for example, that a realistic machine should be able to 
outperform this estimate by over an order of magnitude 35 1. 

A glance at Fig. [7] reveals that it generally adheres to the expected correlation of decreasing scales with increasing 
LNV rates, similar to what is observed for other LNV observable {e.g., Fig. H]). The similarities between the different 
processes extend beyond mere trends to the specific ordering of the operators within each histogram. This reflects 
the common underlying interactions that drive these processes. The operators on the far right of the plot, topping off 
the highest cross-sections, are exactly those operators with the largest m^f , now "split" into three different bars. The 
large bar just below 10^ fb is composed of sub- TeV scale operators with tree-level diagram ZJg-like fermion content. 
Slightly smaller are the expectations for Oig, again dominated by diagram Dg, but characterized by a slightly larger 
Ay scale (around one TeV). Moving down in cross-section, this is followed by the low cutoff scale operators ©54,, ^ and 
©70, dominated by a combination of diagrams Dq and D7. On the opposite end of the plot we point out the large 
bar below lO"^'"^ fb, composed mainly of operators associated to high cutoff scales (Aj^ > 25 TeV). The contributions 
of these operators are dominated by light Majorana neutrino exchange, but their histogram bar contain far fewer 
models than their [3l3Qv counterpart, as many of the latter have been driven up due to new diagram and Z^s 
contributions. In general, the large center-of-mass energies tend to magnify differences between interaction rates that 
were not relevant in low-energy observables. This naively suggests that high energy probes have a higher potential 
for distinguishing different models. 

There are eleven operators that lead to an observably large (as defined earlier) e~e~ — > qqqq cross-section at the 
ILC. Note that all of these were already "ruled out" by current l3(3Qv searches. As discussed in Sec. IIII Al however, 
these bounds only effectively limit the couplings of the new physics to first generation of quarks and leptons, and 
hence, if such a scenario is realized in nature, one should still expect large contributions from decay modes that lead to 
second and third generation final state quarks. In fact, even one such heavy quark is enough to bypass the constraints 
from PPQv for several effective operators. Such reasoning implies that constraints on the new physics flavor structure 
can be made quite strong at a linear collider via analyzes of the flavor of the final state quarks. By identifying and 
comparing the outgoing quark flavor one can extract individual limits on quark-lcpton coupling constants within 
the operators. Additionally, kinematics can be used as a further operator probe. For example, one can potentially 
determine the dominant underlying LNV diagram (say Dq, or Dg) by checking whether the various kinematic 
distributions are characteristic of M^-boson exchange. 

The ILC can cleanly select or discard some LNV scenarios. This characteristic is further enhanced by considerations 
of initial electron polarization. Planned linear colliders have the ability to produce partially polarized beams (80% 
polarization for e~, 40% for [35i l37j). The power of a high energy polarized e~e~ beam is in model identification 
and rejection. Of all operators that yield observably large cross-sections, the e^e^ mode can only probe ©53, and 
therefore any positive LNV signal cleanly identifies this as the operator chosen by nature. In a similar way, the ILC 
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FIG. 7: Distribution of total cross-section for the process e~e~ <l<iq<i and no missing energy at an e~e~ collider with 1 TeV 
of center-of-mass energy. Estimates were obtained assuming the scales derived in Sec. [Ill as well as order one coupling 
constants. The histogram bars are labeled with operator names and color-coded by cutoff scale. Also shown (broken 
vertical line) is the reach of such an experiment assuming 100 fb~^ of integrated luminosity. The vertical axis is truncated to 
best display the relevant features of the plot. 



running in its e^e^j mode can easily observe LNV from O19, C54„, C54dj C'59 and Oeo! ^^nd to a lesser extent, operators 
054(j and O70, and possibly even ©54^. Finally the e^e^^ mode can probe operators Osg, O37 and O^s- Within this 
framework, any LNV detected in one ILC polarization mode will generally not be seen in the others. This statement 
also applies to resonantly enhanced low scale operators that lie outside the observability window. 

While e~e~ collisions only probe effective operators that "talk" to first generation leptons, there are several lepton 
collider processes that allow one to explore other members of the charged lepton family. Future high energy muon 
colliders [s^ could, in principle, also be used to study LNV. In this case, all of the preceding discussions regarding 
the ILC are applicable. Electron linear collider facilities can also be used to study 76" and 77 collisions [39|. 76" 
collisions can be used to probe 76" ^ i'^ + X (and hence the "ea" structure of different LNV operators), while 
77 ^t^a + probes all the different a, f3 charged lepton flavors. For 77 collisions, for example, considering 
projected ILC-like collider parameters, one would expect the same operator distribution as Fig. [7l shifted down in 
cross-section by, roughly, a factor of ^ 10"''. Thus, a handful of operators should be testable at a future 77 collider 
assuming 100 fb"^ of integrated luminosity. 

The preceding analyses carry over to the case of hadron colliders, such as the LHC, in a relatively straightforward 
way. The LHC, or Large Hadron Collider, is a proton-proton machine that will operate at a center-of-mass energy of 
14 TeV and a characteristic integrated luminosity around 100 fb"^ [34] (in its high luminosity mode). The relevant 
LNV variants of Eq. (jlll.ip are dd t~i^uu and uu i'^t^dd with no missing energy. Of course, at center-of-mass 
energies well above a TeV, the proton-proton collisions are dominated by the gluon content of the proton, so most 
interactions at the LHC will be initiated by gluon-gluon and gluon-quark scattering. The dominant LNV subprocesses 
are qg — *■ ia^^QQQ ^^i^ 99 ~^ ^^'^ ^^'^ illustrated in diagrams (a) and (6) of Fig. [51 respectively. These are 

characterized by similar final states as the quark-quark scattering reactions but, given that there is no explicit gauge 
boson field content in the LNV operators in question (Table|I|, their amplitudes are proportional to unimportant order 
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(a) (6) (c) 

FIG. 8: Parton level gluon-gluon and gluon-quark LNV interactions relevant at high energy hadron colliders. Each of these 
yields a same sign dilepton signal with jets and no missing energy. Notice that the final state flavor structure is completely 
arbitrary under the assumption of random order one coupling constants. 



as and coefficients, respectively. The parton level diagram (c) shows the related process gg f.aVp + qq. The rate 
for this process can be estimated, relative to its four jet cousins, by exchanging a final state phase space suppression 
for a single loop suppression. In all three diagrams depicted in Fig. [8l the LNV interaction regions represented by 
large grey dots contain all of the diagrams discussed earlier, meaning that the operator amplitudes calculated for the 
ILC can be recycled in this analysis. While all three bare diagrams are characterized by rates of the same order of 
magnitude diagram (c) leads to missing transverse energy and potentially undetermined final-state lepton number, 
rendering it a less than optimal experimental search mode. Note that, in all of these cases, the external, and internal, 
fermions outside of the LNV interaction region can be of any flavor. Therefore, hadron collider experiments have, 
in principle, access to all LNV operator parameters. Cleanly identifying and constraining all said parameters should 
prove quite difficult for all but the most obvious signatures. The above statements regarding signals at the LHC are 
also applicable at the Tcvatron with some minor, but important, modifications. The Tevatron's pp collisions are at a 
much lower center-of-mass energy, roughly 2 TeV, while the total expected integrated luminosity, less than 10 fb~^ per 
experiment, is orders of magnitude smaller. These factors lead to much lower amplitudes, reduced by approximately a 
factor of {QTcvatmn/ Q LHc) ~ lO"^.-'' The smaller center-of-mass energy also limits the Tevatron's ability to directly 
produce new physics states. With this in mind we conclude that the Tevatron has little or no chance of discovering 
LNV (within this minimal framework). 

A detailed set of predictions for the LHC would require a much more refined analysis, including the effects of parton 
structure functions, flux distributions, and backgrounds, and as such is beyond the scope of this general survey. We 
would, however, like to point out that some of the reactions outlined here are subject to large background rates. While 
SM processes are lepton number conserving, many can fake the LNV signals in the complicated environment of a high 
energy hadronic interaction. The requirement of no missing flnal state energy is particularly hard to accommodate as 
some energy is always lost down the beampipe. As is typically done, one must rely on the less restrictive conservation 
of transverse momentum in order to constrain invisible states, such as neutrinos. SM same-sign dilepton production 
processes arising from, say, M^-boson pair production, are serious potential sources of background. Furthermore, it is 
impossible to predict correlations among final state jets without selecting a particular operator and underlying model 
of new physics, making it difficult to impose general cuts to reduce other hadronic backgrounds. Of course, some of 
the low scale LNV operators yield large enough total cross-sections that even crude analyses may suffice to reveal 
their existence. We conclude by pointing out that a large amount of recent work has been dedicated to LNV searches 
at collider facilities [i^. Most of these approach the subject from the perspective of sub-TeV mass, mostly sterile 
Majorana neutrinos that mix with the active neutrinos and are thus related to light neutrino masses via the seesaw 
mechanism Q. This amounts to one example that leads to the dimension-five operator Oi, but where one assumes 
that the propagating degrees of freedom are twelve or thirteen orders of magnitude lighter than the ultraviolet cutoff 
scale A^} In this case, LNV interactions are dominated by diagram of Fig. [3] (where heavy (weak scale) neutrinos 
are also exchanged), and as such one should make use of specific kinematic cuts to reduce background rates. These 
cuts, however, may also remove LNV signals resulting from many of the scenarios explored here, particularly those 



t Strictly speaking one must also account for the proton's structure functions at the Tevatron's energy scale. Unlike the LHC, where 
collisions are dominated by gluon-gluon interactions, proton collisions at the Tevatron are dominated by valence quark interactions. 
These considerations do not affect our conclusions. 

§ This can be achieved in two different ways. Either the new physics is very weakly coupled, or the new physics - SM couplings are 
finely-tuned [f4| . In order to observe right-handed neutrinos in colliders, the latter must be realized. 
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whose rates are dominated by Dg at tree-level. We urge experimentalists to account for this possibility while analyzing 
future data sets. 



Table H] contains predictions for all the entries rria/j, the Majorana neutrino mass matrix. These are computed 
in the weak basis where the weak interactions and the charged-lepton Yukawa couplings are diagonal, so that the 
eigenvalues of the neutrino mass matrix are the neutrino masses (bounded by oscillation experiments and, say, precision 
measurements of tritium beta-decay [7]), while its eigenvectors determine the neutrino mixing matrix, constrained 
mostly by oscillation experiments. Since different LNV effective operators predict different flavor-structures for the 
neutrino mass matrix, there is the possibility to constrain the different scenarios with existing oscillation data 0]. 
While we can only predict the values of rriafs within, at best, an order of magnitude, it is still possible to extract 
useful information from the derived large scale structure of the expressions. In particular we can test the hypothesis 
of whether A values associated to different lepton flavors are allowed to be of the same order of magnitude. In order 
to obtain more accurate predictions and further probe the fine details of lepton mixing one must succumb to specific 
models, beyond the scope and philosophy of this analysis. 

The mass matrix for the three light Majorana neutrinos can be reconstructed from nine observables: three masses 
mi,m2,TO3, taken to be real and positive; three (real) mixing angles ^12,^23,^13; and three CP-violating phases 
(5, 02, 03- Here, ^ is a so-called Dirac phase that is generally present in the system regardless of the neutrino's nature 
(Majorana or Dirac fermion), while 0i, 4)2 are so-called Majorana phases, only present if the neutrinos are Majorana 
particles (which is the case of all scenarios under consideration here). Oscillation data determine with relatively good 
precision ^12, ^23, Am^2 = "^2 ~ "^1 ^'^^ |A'^i3l = 1"^! ~ '^il- define neutrino masses such that mi < m2 and 
Am^2 < lAmfgl, so that the sign of Amjg remains as an observables which characterizes the neutrino mass hierarchy 
("normal" for Amfg > 0, "inverted" for Amfg < 0). See, for example, [H for details. As for the third mixing angle, 
sin^ 9i3 is constrained to be less than 0.025 (0.058) at 2a (4(t) from a three neutrino global oscillation analysis A 
considerable amount of uncertainty remains. In particular we have only upper bounds on the absolute neutrino mass 
scale, from kinematical measurements such as tritium beta decay [1, Q , plus cosmological observations ^ . Finally, 
the three CP violating phases are completely unconstrained, and we have no information regarding the neutrino mass 
hierarchy. 

The above exp erimental results allow for several different "textures" for ma/i in our weak basis of choice (see, 
for example, [4l|). The purpose of this section is to discuss whether any of the textures predicted by the different 
LNV effective operators is "ruled out" by current observations. Most of the analyzed operators imply "anarchic" Q 
neutrino masses. This simply means that all elements of the neutrino mass matrix are uncorrelated and of the same 
order of magnitude. This hypothesis is known to "fit" the current data very well Q. It will be further challenged by 
searches for 613 (the anarchic hypothesis favors large 6*13 values) and probes that may reveal if the neutrino masses 
are hierarchical or whether two or three of the masses are almost degenerate (anarchy naively predicts the former). If 
future data strongly points towards non-anarchic ruafs, we will be forced to conclude that there is nontrivial "leptonic" 
structure in the dimensionless coefficients A of most of the LNV operators considered here. 

Many of the operators associated with a low neutrino-mass related cutoff scale (A^ < 10 TeV), on the other hand, 
naively predict more structured neutrino mass matrices. Operators 



which radiatively generate neutrino mass elements proportional to distinct charged lepton Yukawa coupling (j/e, 2/^,2/,-), 
yield mass matrices m such that 



Additionally, models described at low energies by O36, 037 and ©38 generate neutrino masses proportional to both 
associated charged Yukawa couplings, such that 



IV. NEUTRINO MIXING 



C7, Cs, Ci9, C20, C'34, C35, C'54„,t,.e,di '^sSa.t,,: , CsG, ^s?, C58, C59, C'eo, O70, C75, 



(IV.l) 




(IV.2) 




(IV.3) 



The strongly hierarchial nature of the charged lepton masses (t/e ^ y^t ^ J/r), implies that the niais elements of 
Eqs. I|IV.2|) and pV.3|) are expected to be hierarchical as well. In particular, the ee matrix element, iriee, proportional 
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to j/e or y1 is, for all practical purposes, negligibly small* in both of these cases. On the other hand, it is well known 
that only a normal neutrino mass hierarchy is consistent with vanishing niee [13 > so that both Eqs. (|IV.2[) and (|IV.3p 
predict the neutrino mass ordering to be normal. In the absence of extra structure, scenarios characterized by the 
LNV operators listed in Eq. pV.ip plus Oag, O37 and Oag will be ruled out if future data favor an inverted mass 
hierarchy, or if the neutrino masses end up quasi-degenerate (regardless of the hierarchy) . As will become clear shortly, 
Eqs. pV.2p and pV.Sp predict that the lightest neutrino mass (mi in this case) is small (< ^JKrn^). 

A more detailed analysis reveals that naive expectations from Eqs. (|IV.2p are already disfavored, while those from 
Eqs. (jlV.SP are virtually excluded. Assuming the normal hierarchy and very small mee, one can find a relation 
between the neutrino mass eigenstates and the oscillation parameters, thus reducing the number of free parameters 
in the mass matrix by one. Consider the diagonalization of the neutrino mass matrix defined by rriap = UM^U'^ 
with = diag(mi, m2e^*'^^ , mae^*"^^) and U the neutrino mixing matrix, expressed in the PDG parameterization. 
In this case, 

mee = mi cos^ 6*12 cos^ 6*13 + m2 sin^ 6*12 cos^ Oize^^'^^ + ^3 sin^ 6'i3e^'("*-'^'^) . (IV.4) 

Setting mee — 0, one can solve for mi and one of the Majorana phases. Recalling that, for the normal mass hierarchy, 
m2 = \/ m\ + Ami2 and raj, = m\ + Am^g, and assuming small ^13 and r/ = Am^2/Am^3, 



mi sin^ 9s „2 '^os 



2 , 



^ "7/.;^. - ^?37^ cos[2(03 - 5)] 



'cosi/22% ^^cos2% 

h^lA arctan (1^^^ sin[2(03 - 5)]) . (IV.5) 
2 2 \ 77 sm 29s J 

One can easily obtain approximate expressions for the other neutrino masses (m2,m3) and hence all elements rriap- 
Upon substituting the numeric best fit oscillation parameters to avoid introducing a needlessly cumbersome expression, 
we get 

= O.Se*^'^^ ( 1 1 1 + 0.716li3e-*('^"2'^-^^ ( 1 1 + 0.4577 ( -1.3 -1 0.61 
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-0.67 

a2 




+ 0.9161^3 cos[2(J- (^3)] ( 1 0.12 0.02 | + 1.276*^3 sin[2(03 - J)] ( 1 1.2 -0.83 (LV.e) 

-0.83 0.56 

Eq. pV.6p suggests a clear hierarchy among the mixing matrix elements. The four, lower box-diagonal /i — t 
elements dominate, followed by the off-diagonal e/i and e/z entries, and finally the vanishingly small mee- Except 
for the vanishingly small mee, which was required a priori^ all of the remaining properties follow directly from the 
experimentally determined mixing parameters. Among the dominant /i — t submatrix, Eq. pV.6p predicts that all 
entries are equal up to small order 77 and 6*13 corrections. The magnitude, and sign, of these "breaking terms" can be 
tuned with the phases 03 and 5, and to a lesser extent by varying 77 and ^13 within their allowed ranges. On the other 
hand, the relative sizes of me^i and mer are expected to be similar but not identical, i.e., nie^ ~ mer ^ {me^ — mer). 

While some of the gross features of Eq. pV.6|l are shared by Eq. pV.2p and Eq. (IIV.3p . a finer analysis reveals 
several disagreements. The major discrepancy lies in the required relations among the matrix elements. Eq. pV.2|) 
predicts that all men- elements are equal, while Eq. pV.3P suggests 77ier *C m^,- ^ mrr- Both of these contradict, 
in different ways, the experimental constraint mer ^ ^/ir ~ m^-T-. Additionally, both Eq. pV.2p and Eq. (IIV.3P 
predict ruee ^ m^^ <C rrirT, while observations require mee ^ "t,^^ « mrr- Similarly, both sets of operators suggest 
niefi mrier while, experimentally, they are constrained to be similar. 

In order to quantify how much Eq. pV.2p and Eq. pV.Sp (dis) agree with our current understanding of neutrino 
masses and lepton mixing, we numerically scanned the allowed mass matrix parameter space assuming the normal 
neutrino mass hierarchy and constraining |mee| < Ue/VT x 1 eV ~ 10~^ eV. It should be noted that, according to this 
relation, niee is allowed to deviate by nearly a factor of ten above naive expectations from mass matrix Eq. pV.2p . thus 
accounting for the possible order of magnitude uncertainties in operator scales and coupling constants. This feature 
is only included for completeness, as one expects that such mee excursions from zero will generally have negligible 



Quantitatively, in the scenarios under investigation, mee values are, respectively, up to order one corrections, yejyr ~ 10 ^ and 
y^/y^ ~ 10~^ times the characteristic mass scale of the mass matrix. 
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FIG. 9: Scatter plots of the symmetric Majorana neutrino mass matrix elements normalized to m-rT- Each panel is produced 
assuming the normal mass hierarchy and parameter constraints insuring that mee < 10"** eV. The light grey region is calculated 
allowing all mixing parameters to vary within their respective 95% confidence intervals. In the purple (darker) regions, the 
solar and atmospheric parameters are held constant while all phases are scanned within their physical ranges and On is varied 
between zero and its 4 a upper bound. The sin'^ On variation is illustrated by varying the shading from dark to light. Also 
indicated by red (closed) and blue (open) dots are the expectations derived from Eqs. (|IV.2p and (|IV.3|) . respectively, along 
with a listing of their associated coordinate values. 



effect on the mass matrix due to the robust nature of Eq. pV.6p . Fig. [9l a scatter plot of mixing matrix elements, 
depicts the result of such a scan. Note that we plot the mass ratios with respect to assumed-to-be-dominant rrirT 
element. The light grey regions of the plot were produced allowing all oscillation parameters to vary within their 
95% confidance bounds and phases to vary within their entire physical range subject to the constraints discussed 
above. In the purple (dark) region, the phases and reactor mixing angle Oi^ are allowed to vary while all other mixing 
arameters are held fixed at their best fit values. We depict the sin^ 9i3 variation from zero to 0.06 (4ct upper bound 
) by varying the purple shading from dark to light. It is easy to check that the numeric (Fig. [S]) and analytic results 
(Eq. (jIV.6P ) are consistent both qualitatively and quantitatively. 

Fig.[9]also depicts the predictions from Eq. (IIV.2P and Eq. pV.Sp with red (closed) and blue (open) dots, respectively. 
As expected, all the predictions from Eq. pV.Sp fall near the origin in each panel and are safely excluded. Because 
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expectations from Eq. (IIV.3P for all map/rriTT are much smaller than one, we also include the dot coordinate values 
for both textures within the figure. In order to render the neutrino mass matrix predicted from O36, O37 and 
©38 consistent with experimental constraints on neutrino masses and lepton mixing, one is required to choose very 
hierarchical A coefficients. In more detail, one needs to choose A values so that all mixing matrix elements are enhanced 
relative to the dominant mrr oc UtDt by numerical factors that range - for different entries - from 100 to 10"". A 
possible mechanism for achieving this is to suppress third generation couplings to new physics, thus driving up the 
ratio map/rrirT along with the required cutoff scale A^. This procedure would have to be accompanied by a more 
modest reduction of the couplings of second generation fermions. Basically, we need to impose a flavor structure 
that "destroys" the naive flavor structure induced by the charged lepton Yukawa coupling hierarchy. We can safely 
conclude that 036, ©37, and ©38, which suggest that the neutrino mass matrix has the form Eq. pV.Sp . are strongly 
disfavored by current neutrino oscillation data and, if somehow realized in nature, must be accompanied by a very 
nontrivial flavor structure. 

On the other hand, the operators listed in Eq. (jlV.ip . which predict Eq. pV.2p . are not quite as disfavored. In this 
case the hierarchies among different mass matrix elements are softer, and one can ask whether the red dots in Fig.[9]can 
move toward the experimentally allowed regions with order 1-10 relative shifts. Many of the predictions are already 
in agreement with experimental constraints, or at least close enough to be easily "nudged" toward acceptable levels 
with order one coefficients. The figure reveals that only m^^ is predicted to be relatively too small. By enhancing 
it by a factor of order yT/y^l ^ 20 one obtains moderately good agreement between Eq. (jIV.2p and experimental 
requirements. We therefore conclude that operators listed in Eq. (jIV.ip are at least marginally allowed by neutrino 
mixing phenomenology. 

While essential for a complete understanding of neutrino masses and mixing, improved measurements of the already 
determined mixing angles and mass-squared differences will not help to further constrain/exclude any of the LNV sce- 
narios in question. Considering our parameter flexibility, only future neutrino experiments that provide qualitatively 
new results can aid in this endeavor. In particular, the experimental determination of the neutrino mass hierarchy is 
essential in order to properly test the scenarios highlighted in this section, as they all predict, in the absence of very 
non-trivial flavor structure in the LNV sector, the normal hierarchy. Next-generation neutrino oscillation experiments 
are expected to provide non-trivial information regardingthe neutrino mass hierarchy. Most rely on a ncutrino/anti- 
neutrino oscillation asymmetry via Earth matter effects [l|, , and depend heavily on a sufficiently large ^13 mi xing 
angle. The possibility that ^13 is vanishingly small, where the standard approach is ineffective, is addressed in (43| 




considering both oscillation and non-oscillation probes. In that case, one can hope to discern the neutrino mass 
spectrum in future neutrino factory / Superbeam experiments coupled with improved constraints on the ef- 
fective masses extracted from tritium beta decay 0] and cosmology Q.^ Note that these non-oscillation probes can 
be independently used to constrain LNV models, as they provide information regarding the the magnitude of the 
lightest mass eigenstate (toi [7713] in the case of normal [inverted] hierarchy). For example, if either cosmological ob- 
servations or tritium beta decay experiments see evidence for non-zero neutrino masses (in more detail, they constrain 
E = and ~ Si'^ilf^eiP respectively) such that T, 0.05 eV or to^^ ^ 0.01 eV, one would conclude, 

assuming a normal mass hierarchy, that mi 3> y^Aro^. This would destroy the possibility of negligibly small mge, 
and hence disfavor the operators that lead to mass matrices of the type Eq. (|IV.2p and (lIV.Sp . Currently, S and 
rriu^ are bounded to be below 0.94 eV and 2.0 eV, respectively, but the sensitivity to these observable is expected to 
significantly improve with next-generation experiments to 0.1 eV [44] and 0.2 eV (45j . respectively. 

V. PHENOMENOLOGICALLY INTERESTING OPERATORS: SAMPLE RENORMALIZABLE MODEL 

Having superficially surveyed a large set of LNV operators, we are now in a position to identify operators with 
"interesting" phenomenological features for further detailed study. One subset of potentially interesting operators is 
characterized by those that, when required to "explain" the observed neutrino masses, are accompanied by a low cutoff 
scale of, say, less than several TeV. Further requiring a small enough rrigf in order to evade current (3/30iy constraints, 
this set contains only seven elements: 017,018,034,035,056,057,058- Of these, all but operators O35 and O58 
(which lead to the zeroth-order neutrino mass matrix Eq. (jIV.2p and a suppressed Wee) should provided a positive 
LNV signal in the next round of double-beta decay experiments, baring specific fiavor symmetries or finely-tuned 
couplings. Furthermore, O56 leads to a /3/30z/ rate that is higher than what is naively dictated by the values of the 



t One traditionally includes the effective /3(30v mass mee given by Eq. I IIII.2II in a neutrino mass hierarchy analysis. However, as discussed 
in Sec. IIII Al m'^ is a potentially convoluted process-dependent quantity that generally has little (directly) to do with neutrino masses. 
For this reason, PfSOu constraints cannot be used to determine the neutrino mass spectrum from the point of this analysis. 




30 



neutrino masses. Finally, with the possible exception of O^e which may mediate observable LNV processes at high 
energy colliders, none of the seven operators above are expect to mediate LNV violating phenomena (as defined here) 
at accessible rates. 

An "orthogonal" subset consists of the higher dimensional operators already "excluded" by l3(3Qv. Not 
including those operators severely constrained by Icpton mixing in Sec. IIVI this list contains 11 elements: 
Ci6, Ci9, ^53, C'54a,b,c,(i, ©59, Oeo, C70, ^Ts- Most of these are associated to cutoff scales of order the weak scale, 
which are likely to already be constrained by different searches for new degrees of freedom with masses around 
100 GeV. Even if those are considered to be excluded, OiQ,Oig,OT5 are "safely" shielded from direct and indirect 
non-LNV searches,* while still mediating potentially observable LNV effects at colliders as long as the new physics 
does not couple, to zeroth order, to first generation quarks (in order to evade the /3/30i/ constraints). 

Regardless of whether these different options for the LNV sector lead to observable LNV phenomena, the low 
extracted cutoff scale of all the operators highlighted above implies that new degrees of freedom should be produced 
and, with a little luck, observed at the LHC or, perhaps, the ILC. Furthermore, the TeV scale has already been 
identified as an interesting candidate scale for new physics for very different reasons, including the dark matter puzzle 
and the gauge hierarchy problem. The fact that, perhaps, the physics responsible for neutrino masses also "lives" at 
the TeV scale is rather appealing. 

In order to study this new physics, as already emphasized earlier, ultraviolet complete manifestations of the physics 
that leads to the effective operators are required. Here we discuss one concrete example. Other examples (for different 
effective operators) were discussed in ^^]. Given a specific LNV operator, it is a simple matter to write down equivalent 
renormalizable Lagrangians. We briefly illustrate this procedure by constructing a renormalizable model that will lead 
to the dimension-eleven operator O^q. It is among the interesting LNV effective operators of the sample highlighted 
above, since it is currently unconstrained by j3(3Qv searches regardless of the quark-flavor structure of the operator, 
while TTigg ^ rriee for PPOv. On the other hand, for O56 is very low (below 500 GeV), so that the new degrees of 
freedom may already be constrained by, for example, Tevatron or LEP data. We will not worry about such constraints 
henceforth, but will only comment on possible phenomenological problems. 

©56 can be accommodated by a wide variety of models, as can be seen from its possible Lorentz structures. In 
terms of scalar/tensor helicity- violating bilinears = 1,it^i/, and vector helicity-conserving bilinears Tc = jfi, these 
are 

O56 = {(LT„Q-')(d=F,„d^)(3^F,i^), (iT,g^)(dT,d^)(d^F,i^), (LT,d=)(QJF,d=)(d^F„F), 
(LT,d^) (g^ F,d^) (dTci^) , (LT„d^) (Q^ Fei^) (dTcd^) , (LTcd^) (g^ F„d^) (d=F,i^) , 

(LTcd^)(g^Fci^)(rfT„d^), (LTei^)(g^F,rf^)(dTerf^), (LT,i^)(g^Fed^)(dT„d^)} x H'^H'e.keji- (V.l) 

It is clear from the chiral field content that these operators depend on combination of helicity-conserving and 
helicity-violating interactions. In particular, it is impossible to form any of the operators in this long list with only 
the addition of vector boson states: new heavy scalar and/or tensor particles are probably required if O5Q is the proper 
tree-level manifestation of the LNV physics at low-energies.^ Furthermore, the couplings of the new physics fields with 
one another must be constrained in order to "block" the presence of lower-dimensional tree-level effective operators. 
This usually implies the existence of new exact (broken) symmetries to forbid (suppress) particular interactions. 

Certain Lorentz structures, those containing only F^, bilinears, can be realized assuming that the LNV ultraviolet 
sector contains only heavy scalar fields and we concentrate, for simplicity, on this possibility (46| . Simple scalar 
interactions that can lead to Ose are shown in the diagram in Fig. [TOl Specifically, these yield the effective operator 
Lorentz structure {UQ^){(P(F){e'^dP)H'' eiktji with the introduction of four charged scalar fields, ^i, (/)2, (/)3, ^4. The 
gauge structure is such that, under (S'C/(3)c, S'C/(2)l, C/(l)y)' , 01 transforms as a (3,3,-1-1/3), (/)2 as (3,1,-2/3), 03 
as (3, 1,-4/3), and 04 as (3, 1,-2/3). While 02 and 04 have identical gauge quantum numbers, they have different 
baryon number (2/3 versus —1/3). 0i has baryon number —1/3, while 03 has baryon number 1/3. Lepton number 
cannot be consistently assigned as it is explicitly violated by two units. 

04, which does not couple to any of the SM fermions, plays an essential role. It acts as a selective "insulator" that 
connects the various interaction terms in such a way as to only alow certain tree-level higher dimensional SM effective 
operators. All renormalizable theories that lead to only very high dimensional effective operators contain one or more 
of these "hidden sector" fields. Note that the new scalar fields should not acquire vacuum expectation values in order 



* Generic new degrees of freedom at the weak scale are constrained by direct and indirect searches at high energy colliders (e.g., resonances 
and effective four-fermion interactions, respectively), flavor-violating {e.g., fi e-y), and high precision experiments {e.g., measurements 
of the anomalous muon magnetic moment). 

t Other possibilities include heavy vector-like fermions. 

•t In the case of U{1}y, 'transfoms as X' means 'has hypercharge X'. 
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FIG. 10: Sample scalar interactions that lead to the "interesting" effective operator Oz,6 with the Lorentz structure 



to avoid the presence of lower dimensional irrelevant operators that are likely to dominate low-energy phenomenology 
and - much more important - to prevent the spontaneous breaking of color or electromagnetic charge. 

Given the scalar field content as well as its transformation properties under SM global and local symmetries, it is a 
simple matter to write down the minimal interaction Lagrangian density for the system. A candidate renormalizable 
Lagrangian is 

C = + {\DM'' + M^\4>^\'') + yiQL^i + V2d'd^4>2 + Vze'<f4>3 + XuMiHH + A234Af 0203^4 + h.c. . (V.2) 

i 

Each term in Eq. (jV.2[) . including those involving covariant derivatives D^, is implicitly assumed to respect the gauge 
representations of the associated (pi fields, as defined above. The Yukawa-type couplings yi, as well as the Xi scalar 
vertices are dimensionless, and assumed to be of order one, while we assume all scalar masses Mi to be of the same 
order of magnitude. In this case, A ~ Ali. In the A234 term, an overall mass scale M has been "factored out" and 
is assumed to be of the same order as the Mi. Note that we neglect generation indicies, which are implied. In the 
case A234 = 0, lepton number is a classical global symmetry of Eq. (|V.2[) . and one can view this three-scalar coupling 
as the source of lepton number violation. One may even envision a scenario where lepton number is spontaneously 
broken by the vacuum expectation value of some SM singlet ^5 scalar field, ((^5) = M. 

Provided all Mi are around 0.5 TeV, as required if this Lagrangian is to "explain" the observed light neutrino masses, 
LNV is certainly not the only (or even the main) consequence of this model. The yi and ys terms for example, will 
mediate /i — > e-conversion in nuclei at very dangerous levels if their flavor structure is generic. 4)2 can be resonantly 
produced in dd-collisions, while 4>i and ( ^3 q ualify as scalar lepto-quarks, which are constrained by high energy collider 
experiments, including those at HERA [431, to weigh more than a few hundred GeV [l^. For more details, we refer 
readers to, for example, the Particle Data Book TT*! and references therein. 

We will conclude this discussion by adding that several other effective operator can be realized in a very similar 
way. Oi9, for example, if it manifests itself with the Lorentz structure {UQ^){df'df-){e'^u'^)eij , can be realized by a 
Lagrangian very similar to Eq. (|V.2[) where the df^ field in the ?y3-coupling interaction is replaced by a u'^ field, and 
the (pi field is replaced by an SU{2)l singlet (it is a triplet in Eq. (|V.2p ). Of course, hypercharge assignments for the 
(pi also need to be modified in a straight forward way. The associated non LNV phenomenology is similar, except for 
the fact that A^, for Oig (around 1 TeV) is larger than the one for Osg and hence O19 is less constrained by current 
experimental data. On the other hand, Oig predicts potentially much larger rates for LNV observables at colliders 
(see Fig,[71). 

Our definition of "interesting" is arbitrary and motivated only by the fact that the physics of the "interesting" 
operators highlighted earlier in this section will probably be explored at next-generation collider and high-precision 
experiments. One may argue that many operators which lead to the observable neutrino masses for high values of 
Ay are interesting on their own right, either due to the theoretically pleasing properties of their associated potential 
ultraviolet completions, or by some observational peculiarity. There are many examples of the first type ranging 
from the different manifestations of the seesaw mechanism [l^, [ll|, [l2| to the Zee model [i^ and the minimal 
supersymmetric SM with R-parity violation [49] . Dedicated analysis of these cases have been widely pursued in the 
literature and will not be discussed here. We would also like to point out that some effective operators, like O7 and 
Og, are, according to our criteria, very "uninteresting." Both Oy and Og predict unobservably suppressed j3[3Qv rates 
(both predict small mee) and equally hopeless collider prospects given that they are associated to very high cutoff 
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scales, « 4 X 10^ TeV and Aj^ « 6 x 10^ TeV, respectively. If either of these operators are responsible for the 
observed tiny neutrino masses, it is quite possible that we may never directly detect LNV. It is curious to consider 
possible means of indirect detection or other observable consequences of the different ultraviolet completions of such 
scenarios.^ It would also be interesting to ask whether either of these elusive models has any underlying theoretical 
motivation or whether they allow one to solve other outstanding problems in particle physics. 

VI. DISCUSSION AND CONCLUDING REMARKS 

If neutrino masses are a consequence of lepton- number violating physics at a very high energy scale (higher than the 
scale of electroweak symmetry breaking), new physics effects ~ including the generation of neutrino Majorana masses 
- at low enough energies can be parameterized in terms of irrelevant operators whose coefficients are suppressed by 
inverse powers of an effective cutoff scale A. As discussed before, A is, roughly, the energy scale above which new 
degrees of freedom must be observed if the new ultraviolet physics is perturbative (if the new physics is very weakly 
coupled, the masses of the new degrees of freedom can be much smaller than A). We have explored a very large class 
of such scenarios through 129 irrelevant operators of energy dimension less than or equal to eleven that violate lepton 
number by 2 units. These are tabulated in first two columns of Table U along with a summary of our results. 

Analyzing each effective operator individually, we estimated the predicted general form of the Majorana neutrino 
mass matrix. Our results are listed in the third column of Table HI By comparing each such estimate with our 
current understanding of neutrino masses, we extracted the cutoff scale h-u of each effective operator, assuming that it 
provides the dominant contribution to the observed neutrino masses. These results are listed in the fourth column of 
Table [J assuming light neutrino masses equal to 0.05 eV (the square root of the atmospheric mass-squared difference), 
and are summarized as follows. Depending on the field content and dimension of the irrelevant operator, the "lepton 
number breaking scale" A^ is predicted to be anywhere from the weak scale (~ 0.1 TeV) all the way up to 10^^ TeV 
(see Fig. [2]). This means that, depending on how lepton number is violated and communicated to the SM, the mass 
of the associated new degrees of freedom is predicted to be anywhere between 100 GeV and 10^^ TeV, even if all new 
physics couplings are order one. We note that in the case of all variations of the seesaw mechanism (Oi), neutrino 
physics constrains = 10^^ TeV such that the new degrees of freedom are either unobservably heavy, extremely 
weakly coupled, or their couplings to the SM degrees of freedom are finely-tuned. It is fair to say that this behavior is 
not characteristic of all LNV ultraviolet physics. One sample ultraviolet theory that leads to dimension-eleven LNV 
effective operators was discussed in Sec. [V] Other examples can be found in [2|, and include supersymmetry with 
trilinear R-parity violation and the Zee model. 

Assuming that a particular operator is responsible for nonzero neutrino masses, it is straight forward to ask whether 
it leads to other observable consequences. Here, we concentrated on several LNV observables, and included future 
LNV searches at the LHC and future lepton machines (like the ILC), along with their ability to directly produce 
(and hopefully observe) new physics states lighter than several hundred GeV. In column five of Table [H we list the 
most favorable modes of experimental observation for each operator. The different relevant probes are: neutrinoless 
double-beta decay (/3/30i^), neutrino oscillation and mixing (mix), direct searches for new particles at the LHC (LHC) 
and ILC (ILC), and virtual LNV effects at collider facilities (HElnv). We find it unlikely that other probes of LNV, 
including rare meson decays, should yield a positive signal in the forseeable future. This conclusion is strongly based 
on the fact that, for all of our analysis, we assume that all new physics degrees of freedom are heavier than the weak 
scale. While the vast majority of operators is most sensitive to searches for neutrinoless double-beta decay, that is not 
true of all operators. Some lead to relatively suppressed rates for PfiOv (mostly because they lead to mass matrices 
with a very small niee) even if they are associated to Ai, < 1 TeV, indicating that, for these scenarios, we are more 
likely to observe the physics behind neutrino masses directly at colliders than to see a finite lifetime for l3j3Qv. Other 
scenarios naively lead to l3l3Qv rates orders of magnitude higher than what is currently allowed by data. If these are 
responsible for the generation of neutrino masses, the new physics is constrained to be somewhat decoupled from 
first generation quarks (for example). In this case, there is hope that LNV phenomena at colliders, which are not 
restricted to first generation quarks, occur with non-negligible rates. 

The sixth column of Table [J lists the current "status" of the operator as either experimentally unconstrained 
(U), constrained (C), or disfavored (D). Such labels are assigned based only on the experimental probes reviewed 
in this work. By arbitrary convention, an 'unconstrained' operator can safely accommodate all existing data even 
if one assumes all its flavor-dependent dimensionless coefficients to be of order one. A 'constrained' operator can 



S This is very similar to the case of d. The main redeeming feature of Oi, other than its simphcity, is the fact that many of its ultraviolet 
completions allow one to explain the matter— antimatter asymmetry of the universe [T3| . 



33 



accommodate all existing data after one allows some of the different flavor-dependent dimensionless coefficients to be 
suppressed with respect to the dominant ones by a factor of 100 or so (as described above). 'Disfavored' operators can 
only accommodate all data only if "tuned" much more severely than the 'constrained' ones, and are usually in trouble 
with more than one "type" of constraint. A glance at column six reveals that 11 out of the 129 operators are disfavored 
by current data. The most stringent constraints come from [3(30^, while all 'disfavored' operators are associated to 
cutoffs at or below 1 TeV. Three of the 'disfavored' operators, Oae, O37, and Oas, are also in disagreement with the 
neutrino oscillation data (see Sec. II V[) . 

Our results illustrate that, as far as "explaining" neutrino masses, the model-building scene is wide open even if 
one postulates that neutrino masses arise as a consequence of lepton number violating, "heavy" physics. Significant 
progress will only be achieved once more experimental information becomes available. The observation that neutrino- 
less double-beta decay occurs with a nonzero rate will help point us in the right direction, but will certainly not reveal 
much about the mechanism behind neutrino masses. A more complete picture can only arise from combined infor- 
mation from several observables, including other LNV observables and the search for new physics at the electroweak 
scale. Other important experimental searches, not discussed here, include all lepton-number conserving "leptonic" 
probes, such as precision measurements of the anomalous magnetic moment of the muon, searches for leptonic electric 
dipole moments, searches for charged-lepton flavor violation, and precision measurements of neutrino-nucleon and 
neutrino-lepton scattering. 
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